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ABSTRACT 


I 


Experiments  and  analyses  at  die  Applied  Physics  Laboratory  and 
elsewhere  have  shown  that  burning  b  neath  external  surfaces  of  a  super¬ 
sonic  vehicle  can  produce  lift  and  rc-cuce  drag.  One  great  advantage  of  ex¬ 
ternal  burning  and/or  burned  gas  expansion  is  that  of  cooling  a  greater  por¬ 
tion  of  the  engine  surfaces  by  radiation  to  the  atmosphere.  This  reduction 
in  internal  cooling  requirements  should  lead  to  a  weight  reduction.  A  corol¬ 
lary  advantage  in  increased  lift- drag  ratio  is  obtained  if  the  burned  gas 
expansion  occurs  beneath  a  rearward- facing  surface  of  a  two-dimensional 
or  half-round  engine. 

Preliminary  analyses  of  some  two-dimensional  external  burning  ram 
jet  (ERJ)  and  external  expansion  ramjet  (EERJ)  models  yield  some  en¬ 
couraging  results.  An  analysis  suggests  that  neat  addition  beneath  the  rear¬ 
ward-facing  surface  of  a  simple  inverted  1  iangilar  section  (ERJ)  in  the 
form  of  a  normal-plane  flame  zone  stabilize  1  oe.iind  a  normal  shock  wave 
would  produce  net  thrust  only  when  the  height  of  the  flame  zone  exceeds  0.  3 
or  0. 4  of  the  expansion  ramp  length.  Such  a  mod'l  is  of  doubtful  practi¬ 
cality,  however,  and  does  not  appear  to  be  representative  of  burning 
patterns  observed  at  APL  and  NASA.  A  model  if  greater  promise  is 
a  compromise  between  a  conventional  ramjet  and  ax.  all- external  ramjet 
which  is  sti)1  designed  to  permit  considerable  heat  rei:ef  for  the  struc¬ 
ture  by  radiation  to  the  atmosphere,  to  provide  lift,  and  to  provide  maxi¬ 
mum  packaging  or  fuel  storage  volume  within  the  engine  body.  Literally, 
it  is  a  two-dimensional,  external  expansion  ramjet  (EERJ)  with  a  short 
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subsonic  duct  (but  very  little  subsonic  diffusion)  and  with  critical  heat 
addition  at  the  duct  exit  serving  as  the  nozzle  "throat.  " 

For  the  EERJ  with  an  ideal,  nonreflecting  expansion  surface, 
maximum  thrust  coefficients  of  0.26  to  0.  90  are  calculated  for  equiva¬ 
lence  ratios  from  0.25  to  1.0  using  kerosene  fuel  at  Mach  5.  Fuel 
specific  impulses  for  high-thrust,  low-lift  geometries  compare  with 
those  calculated  for  a  conventional  ramjet  using  similar  assumptions 
on  diffuser  total  pressure  recovery  (38. 3  per  cent),  combustion  effi¬ 
ciency  (100  per  cent),  and  nozzle  velocity  coefficient  (0.  97  per  cent). 

By  varying  the  cowl-body  geometry,  a  significant  lift  force  can  be  ob¬ 
tained  at  the  expense  of  some  thrust  loss.  It  appears  that  the  best 
overall  performance  will  be  obtained  with  a  flat-top  engine  at  zero 
angle  of  attack  with  the  cowl  and  effusor  geometry  designed  to  turn 
the  burned  gas  back  parallel  to  the  flat  top  while  expanding  it  to  a  pres¬ 
sure  of  approximately  three  times  the  ambient  pressure.  Expansion  to 
lower  pressures  appears  to  be  neither  practical  nor  desirable.  In  order 
to  treat  a  simple  case  of  thrust  equals  drag  for  maximum-range  cruise, 
the  product  of  fuel  specific  impulse  and  overall  lift-drag  ratio  is  esti¬ 
mated  for  the  engine  plus  added  wings  for  which  L ID  =  6.  A  maximum 
Ij.  L/D  overall  of  9000  seconds  is  estimated  for  a  kerosene-air  equiva¬ 
lence  ratio  of  0.  6  at  Mach  5.  It  is  estimated  that  the  conventional  ramjet 
plus  added  wings  (wings  L/D  =  6)  would  have  an  1^  L/D  approximately  8 
percent  smaller. 

Effects  of  vectoring  the  EERJ  nozzle  and  advantages  and  dis¬ 
advantages  of  substituting  a  simple  ramp  in  place  of  the  non-reflecting 
surface  are  explored  briefly.  The  lift  and  range  performance  of  a 
simple  ramp  appears  to  be  always  inferior  to  the  maximum  values 
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for  an  ideal  effusor  because  of  appreciable  refractions  of  the  charac¬ 
teristics  for  the  Prandtl-Meyer  expansion  of  a  real  gas  along  a  flat 
surface. 

A  comparison  cf  the  range  versus  weight  performance  of  mis¬ 
siles  using  ramjets,  ERJ,  or  various  rockets  for  second  stage  propul¬ 
sion  demonstrates  the  great  superiority  of  the  air-breathers,  and 
particularly  of  an  advanced,  well -integrated  ERJ,  in  this  regard. 

A  first  qualitative  consideration  of  design  philosophy  for 
hypersonic  transport  vehicles  incorporating  ERJ's  for  cruise  power 
suggests  that  there  is  apt  to  be  a  variety  of  ways  in  which  the  ERJ 
can  be  integrated  with  suitable  airframes  and  auxiliary  lower  speed 
engines  to  provide  a  practical  hypersonic  aircraft.  Three  possible 
paths  for  the  design  are  pictured  and  discussed  in  general  terms. 
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The  conventional  ramjet  engines  for  the  Mach  4  to  5  speed  range 
of  these  missiles  [see  Fig.  1(a)]  have  just  about  reached  the  limits  of 
structural  design  with  the  more  commonly  available  materials  and  simple 
engine  cooling  techniques  such  as  cooling  shrouds  or  thin  interior 
ceramic  coatings.  An  engine  of  similar  type,  that  is,  with  subsonic  in¬ 
ternal  burning  and  internal  exhaust  nozzle,  designed  for  appreciably 
higher  flight  speeds,  Mach  7  to  10  [Fig.  1(b)]  would  have  enormously 
greater  cooling  problems.  Although  the  combustion  chamber  could  be 
much  smaller  (because  the  heat  release  rate  would  be  much  larger), 
the  subsonic  diffuser  would  still  present  a  large  heat  transfer  area  and 
the  exhaust  nozzle  length  and  surface  area  (supersonic  portion)  would 
be  much  greater  because  of  the  larger  expansion  ratio  required.  It 
has  been  estimated  (Ref.  I)1  that  for  a  Mach  8  engine  the  cooling  re¬ 
quirement  imposed  by  the  exhaust  nozzle  would  represent  at  least  half 
of  the  total  cooling  load  for  the  engine,  and  the  small  combustion  chamber 
would  represent  another  1 5  oercent  of  the  total. 

An  obvious  way  to  reduce  the  amount  of  cooling  required  for  such 
a  hypers?  nic  engine  by  a  factor  near  two,  then,  would  be  to  use  a  plug- 
type,  external  expansion  nozzle.  Moreover,  if  a  short  combustion  zone 
could  be  stabilized  at  the  diffuser  exit  to  thermally  choke  the  exit  and 
serve  as  the  nozzle  ’’throat,  ”  the  cooling  problem  would  be  further 
* 

This  report  was  originally  issued  as  Confidential  CM-948  in  1959.  It 
has  been  considered  of  enough  general  interest  to  warrant  revision  to 
make  it  Unclassified  and  to  reissue  it  as  TG-892. 

References  may  be  found  on  pp  85-88. 
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reduced,  because  the  flame  gases  and  adjacent  surfaces  could  "see"  a 
larger  sector  of  space  and  of  external  reradiating  nozzle.  Such  a  design 
[Fig.  1(c)]  has  been  called  an  "external  expansion  burning  ramjet"  or 
"EERJ."  Of  course,  the  ultimate  design  for  making  maximum  use  of 
cooling  by  radiation  to  the  atmosphere  would  be  the  pure  external  burn¬ 
ing  ramjet  or  ERJ  with  no  ducted  flow  at  all  [Fig.  1(d)],  but,  as  will 
be  shown  later,  it  appears  much  more  doubtful  that  the  latter  would  be 
practical  as  a  primary  power  plant  with  acceleration  capability. 

Another  approach  to  the  hypersonic  ramjet  engine  is,  of  course, 
the  ducted  but  supersonic  burning  concept  [see  Fig.  1(e)].  The  main 
reason  for  wanting  to  keep  the  flow  supersonic  throughout  is  to  avoid  tak¬ 
ing  a  large  shock  loss  in  the  diffuser  by  decelerating  the  incoming  air 
only  to  some  intermediate  super  sonic  Mach  number,  say  Mach  2  or  3. 
Analyses  (Refs.  2  and  3)  suggest  that  with  storable  fuels  (e. g.,kerosene) 
flight  Mach  numbers  in  the  range  from  approximately  5.  6  to  10  or  even 
higher  should  be  feasible  insofar  as  calculated  net  thrust  coefficients 
are  concerned.  Such  a  system  requires  the  initiation  and  completion 
of  combustion  in  a  wholly  supersonic  gas  stream  in  which  the  velocity 
would  be  of  the  order  of  4000  to  9000  ft /sec.  This  means  that  fuel 
atomization,  mixing,  and  combustion  should  be  completed  in  a  few 
tenths  of  a  millisecond.  Reported  experimental  work  on  supersonic 
burning  of  storable  fuels  has  been  limited  to  date  (Refs.  4  and  5). 

The  possibilities  of  burning  behind  a  shock  wave  ("standing  wave"  or 
"detonative"  ramjet)  are  also  being  studied  (Refs.  6  through  10) 
but  there  have  been  no  reports  of  engine-type  experimental  work  in 
this  area.  Further  discussion  of  these  types  is  considered  beyond 
the  scope  of  this  paper.  Another  potential  advantage  of  the  EERJ 
(relative  to  a  conventional  ramjet)  arises  if  one  considers  the 
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simultaneous  production  of  thrust  and  lift.  That  is,  the  expansion  of 
burned  gases  beneath  a  rearward -facing  surface  produces  a  force  on 
that  surface  which  may  be  resolved  into  thrust  and  lift  component 
forces.  Proper  integration  of  such  an  engine  into  a  composite  aircraft 
design  might  therefore  result  in  a  significant  increase  in  lift /drag  ratio 
for  the  vehicle. 

Experimental  work  at  the  Lewis  Research  Center,  NASA, 
produced  some  very  interesting  results  (Refs.  11  through  14).  For 
example,  it  was  demonstrated  that  burning  aluminum  borohydride  be¬ 
neath  flat  plates  and  airfoils  in  a  Mach  2.  5  stream  can  produce  signif¬ 
icant  lift  forces  (enough  to  double  L/D)  and  that  burning  around  an 
ogive- cylinder  body  could  reduce  drag  to  practically  zero.  Linearized 
flow  analyses  include  those  of  Gazley  (Ref.  15),  Mager  (Ref.  !6),  and 
Willmarth  (Ref.  17),  all  of  whom  conclude  that  net  thrust  should  be 
possible  as  a  result  of  heat  additions  to  airfoils  in  supersonic  flight. 
Willmarth  also  considers  the  flow  field  produced  by  large  rates  of  heat 
addition  near  the  under  surface  of  a  flat  plate  at  M  =  2  and  concludes 
that,  although  for  large  heating  rates  (flow  deflections  of~  20°)  the  ef¬ 
ficiency  drops  tc  approximately  half  that  for  the  linearized  analysis,  it 
is  still  a  respectable  efficiency. 

To  the  best  of  the  authors'  knowledge,  none  of  the  above- 
mentioned  investigators  conclusively  demonstrated  by  experiments  that 
positive  net  thrust  could  be  produced  by  external  burning.  The  first 
successful  tests  in  the  experimental  program  at  APL  were  made  in 
March  1959  Aluminum  triethyl  fuel  was  burned  beneath  a  simple,  all¬ 
external,  flattop,  10°  by  10°  wedge  model  in  a  6-  by  7-inch  Mach  5 
tunnel  with  a  stagnation  temperature  of  approximately  1200°F  and  a 
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2 

stagnation  pressure  of  78  lb /in  abs.  Net  thrusts  were  indicated  at 
small  negative  angles  of  attack,  with  the  particular  model  and,  of  course, 
significant  lift  forces  were  produced. 

The  present  paper  is  confined  to  the  preliminary  analytical 
studies  which  have  been  made  at  APL,  together  with  some  discussion 
(contributed  by  R.  H.  Cramer)  of  the  integration  of  two-dimensional 
ERJ  engines  into  hypersonic  airframes. 
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I 

T  II.  THE  ALL-EXTERNAL  MODEL  OF  SMITH  AND  DAVIS 


This  section  presents  a  discussion  of  results  obtainable  from  an 
all- external  model  such  as  proposed  by  Smith  and  Davis  (Ref.  18),  who 
reasoned  that,  although  the  earlier  discussions  of  Pinkel  et  al. ,  (Ref. 

19),  Wald  (Ref.  20),  and  others  regarding  heat  additions  to  supersonic 
streams  were  interesting,  no  practical  way  of  adding  heat  to  supersonic 
streams  was  known.  They  therefore  considered  the  case  where  heat 
addition  would  result  in  breakdown  of  the  supersonic  stream  to  produce 
a  normal  shock.  The  actual  combustion  could  then  be  carried  out  under 
subsonic  flow  conditions.  Their  basic  model  consisted  of  a  normal 
shock  wave  supported  by  a  plane  flame  front  (probably  stabilized  by  a 
grid  flameholder)  which  was  perpendicular  to  a  flat  surface  of  an  air¬ 
foil. 

In  the  present  report,  this  plane  flame  and  shock  system  will  be 
considered  in  relation  to  the  rearward-facing  surface  of  a  two-dimensional 
airfoil  having  the  shape  of  an  inverted  isosceles  triangle  with  forward  and 
trailing  wedge  angle  a.  The  model  is  sketched  in  Fig.  2.  Undisturbed 
air  approaches  through  region  0,  passes  through  the  oblique  bow  shock 
into  region  1,  undergoes  an  isentropic  Prandtl-Meyer  expansion  around 
the  hip  into  region  2,  then  passes  through  the  normal  shock  and  flame 
and  expands  in  region  5.  Station  3  refers  to  stream  properties  imme¬ 
diately  downstream  of  the  normal  shock  and  Station  4  to  properties  just 
upstream  of  the  flame.  The  normal  shock  wave  is  assumed  to  remain 
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Fig.  2  MODEL  OF  ALL-EXTERNAL  RAMJET  WITH  NORMAL  PLANE  FLAME 
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parallel  to  the  flame  front  until  it  intersects  the  bounding  stream  line, 
where  it  abruptly  changes  to  an  oblique  shock  at  the  edge  of  the  line 
flame. 


A  linearized  solution  was  obtained  by  Smith  and  Davis  for  the 
normal  force  on  the  surface  between  the  shock  and  the  flame  which  As 
called  Fgn  herein.  (A  more  sophisticated  treatment  of  the  flow  in  region 
3  and  4  has  just  been  done  by  Woolard  (Ref.  21 )  to  be  discussed  later. ) 
For  region  5,  the  pressure  in  the  burned  gas  was  determined  by  mass 
and  momentum  balances.  Immediately  downstream  of  the  flame,  the 
pressure  in  the  adjacent  air  stream  is  increased  to  by  the  oblique 
shock  which  results  from  its  flow  deflection  through  the  angle  6  by 
the  slip  plane.  The  burned  gas  is  assumed  to  leave  the  flame  at  sonic 
velocity  and  to  expand  supersonically  according  to  the  one-dimensional 
flow  approximation  (pressure  constant  along  any  plane  perpendicular  to 
the  wedge  surface).  The  nondimensionalized  normal  forces  F^/Pgh 
F,.  /P_h  (integrated  from  excess  pressure-distance  curves)  were  given 

OH  u 

by  Smith  and  Davis  for  lean  and  stoichiometric  hydrogen  flames  for 
M^'s  of  2,  3,  and  4,  In  order  to  find  the  minimum  ratio  of  flame  height 
to  expansion  ramp  length,  h/i,  required  to  produce  a  positive  net 
thrust  on  the  airfoil,  these  normal  forces  have  been  substituted  into 
the  following  simple  equation  (see  Fig.  2). 


F  ,  =  P^i  sin  Of 
net  O 


(F  +  F  )  P 
v  3n  5n;  2 


P  -  P 
1  2 


V  ' 


In  Fig.  3,  (h/i)^.^  has  been  plotted  against  flight  Mach  number 
Mq  (where  Mq  <  M^)  for  their  stoichiometric  hydrogen  forces.  The 
wedge  angle  Of  has  been  taken  as  10°.  It  appears  that  h li  should  exceed 
3/4  even  for  flight  Mach  numbers  in  the  range  3  to  5. 


_  q  _ 
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Fig.  3  MINIMUM.  RATIO  OF  PLANE  FLAME  HEIGHT  TO  EXPANSION 
RAMP  LENGTH  REQUIRED  FOR  NET  THRUST 
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Since  the  above  result  does  not  look  very  practical,  an  effort  was 
made  to  determine  the  effect  of  refining  the  calculation  by  applying  the 
method  of  characteristics  to  the  flow  behind  the  flame.  For  simplicity, 
the  "lean  limit"  or  "critical  heat  addition"  case  was  studied,  in  which 
the  separation  between  shock  front  and  flame  is  infinitesimal  and  just 
enough  heat  is  added  to  produce  sonic  velocity  at  the  flame  exit  plane. 
The  requirement  that  =  1  is  the  Chapman- Jouguet  condition,  so 
that  identification  with  the  well-known  model  of  a  detonation  wave  is 
complete  for  this  lean  limit  case.  Pressure  and  direction  of  burned-gas 
and  adjacent-air  flows  were  equalized  by  turning  the  burned  gas  in  a 
Prandtl-Meyer  expansion,  letting  the  expansion  waves  radiate  from  the 
outboard  tip  of  the  flame  to  the  wedge.  It  was  found  that  after  a  few 
degrees  of  expansion,  the  refraction  of  these  Mach  waves  by  the  waves 
reflected  from  the  surface  were  slight,  so  that  refraction  effects  were 
neglected  after  x/h  =  1. 

The  resultant  pressure  distribution  on  the  surface  behind  the 
flame  is  shown  in  Fig.  4  for  the  case  y  =  1.4,  M2  =  6.  03.  This  M2 
corresponds  to  a  flight  of  5  for  a  wedge  angle  of  10°.  A  one- 
dimer.sional  (channel)  flow  solution  for  the  same  conditions  is  included 
for  comparison.  Substitution  of  the  forces  determined  by  the  areas 
under  these  two  curves  into  Eq.  (1)  shows  that  the  characteristic 
solution  gives  a  minimum  h/i  of  0.42  compared  to  1.4  for  the  channel 
flow  solution.  While  the  former  value  is  much  more  reasonable  than 
the  latter,  one  still  wonders  if  such  a  flame  height  could  be  practical. 

Woolard  (Ref.  21)  has  given  a  more  realistic  treatment  to  this 
general  model,  particularly  for  supercritical  h^at  addition  in  which 
the  normal  shock  is  detached  from  the  plane  flame.  Figure  5  is  a 
sketch  of  the  flow  pattern  with  detached  shock,which  shows  that  a 
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Fig.  *1  PRESSURE  DISTRIBUTION  ON  SURFACE  BEHIND  NORMAL  PLANE 
FLAME  CALCULATED  BY  TWO  METHODS 
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Fig.  5  FLOW  CONFIGURATION  FOR  HEAT  ADDITION  GREATER  THAN 
CRITICAL  WITH  NORMAL  PLANE  FLAME  (AFTER  WOOLARD) 
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second  shock  must  emanate  from  the  outboard  edge  of  the  heat  addition 
region  to  se~ve  the  purpose  of  equalizing  flow  direction  and  pressure  be¬ 
tween  the  spilled  air  and  the  burned  gases. 

For  the  region  between  the  detached  shock  and  the  flame,  Woolard 

calculates  normal  surface  force  coefficients  approximately  four  times  as 

great  as  those  given  by  Smith  and  Davis,  and  total  after  forces,  (F_  + 

in 

Fg^l/P^h,  approximately  twice  as  great.  The  figures  tabulated  are  based 


Spillage 


F 


3n 


P2h 


Smith  and  Davis 


3  0.266  0.52  0.39  5.0  3.7  1.1  3.1 

4  0.321  0.42  0,23  4.8  8.8  1.  1  7.2 


on  Smith  and  Davis'  value  for  M4  for  stoichiometric  hydrogen  flames. 

The  h/i  values  correspond  to  those  required  to  bring  Woolard's  C,.  to 

on 

90  percent  of  its  maximum  (expansion  to  i  =  ®)  value;  for  these  h/i's  a  net 
thrust  would  be  obtained  on  a  10°  by  10°  wedge  with  Woolard's  values,  but 
not  with  Smith  and  Davis*  (refer  to  Fig.  3). 

The  results  of  Woolard  therefore  make  the  pure  external  ramjet, 
as  represented  by  this  normal  plane  flame  model,  a  little  easier  to  be¬ 
lieve,  though  the  minimum  h/i  for  net  thrust  would  probably  still  be  of 
the  order  of  0.3  (no  exact  correlation  between  Woolard's  results  and 
those  previously  quoted  has  been  attempted  as  of  this  writing).  Some 
fuel  distribution  system  and  flameholder  extending  into  the  stream  would 
be  required  to  produce  and  stabilize  such  a  plane  flame;  these  components 
would  add  drag  and  increase  the  minimum  h/i.  No  way  has  been 
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visualized  to  circum/ent  the  loss  caused  by  the  Prandtl- Meyer  expansion 
around  the  hip  or  mid-chord  of  an  airfoil- -if  the  flame  is  on  the  rear 
surface,  then  M2  >  M^,  T 2  <  Tq,  P2  <  P^;  is  the  flame  is  located  at  or 
ahead  of  the  mid-chord,  the  forward  normal  force  now  produces  drag 
and  the  burned  gases  take  the  expansion  loss  so  that  the  after  force  (be¬ 
hind  the  flame)  is  largely  dissipated.  Thus,  one  concludes  that  a 
reasonable  compromise  might  be  to  return  to  the  use  of  a  short  sub¬ 
sonic  duct  to  capture  compressed  air  and  turn  it  around  the  hip  before 
adding  heat  and  expanding  the  burned  gases  externally. 
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III.  MODIFIED  MODEL  WITH  SHORT  COWL  (SUBSONIC  DUCT) 
TO  TURN  COMPRESSED  AIR 


As  noted  in  the  previous  section,  a  serious  shortcoming  of  an 
all-external  model  as  a  thrust-producer  is  that  in  order  to  have  the 
burned  gas  expansion  occur  beneath  a  rearward  facing  surface,  the 
flow  must  be  turned  around  a  corner  or  hump  in  the  model.  If  the 
flame  and  its  associated  shock  system  are  located  on  the  rear  surface, 
this  turning  requirement  causes  the  incoming  air  to  be  accelerated  as 
it  goes  through  a  Prandtl-Meyer  expansion  around  the  corner,  so  that 
the  Mach  number  M  approaching  the  flame  shock  is  actually  higher 
than  the  free  stream  or  flight  Mach  number.  This  high  M2  results 
in  a  high  total  pressure  loss  through  the  flame  shock  system 
(P  /P  =  0.03  for  Mach  5  flight  of  the  10c  by  10°  model).  Conversely, 

to  tb 

if  the  flame  is  established  just  upstream  of  the  knee  in  the  model,  then 
that  portion  of  the  burned  gas  expansion  which  corresponds  to  the  turn¬ 
ing  angle  around  the  knee  cannot  exert  any  pressure  on  the  rearward  - 
facing  surface;  moreover  any  excess  pressure  on  the  forward- facing 
surface  caused  by  combustion  contributes  to  drag  rather  than  thrust. 

The  external  expansion  ramjet  avoids  the  flow  turning  loss.  A 
diagrammatic  sketch  of  the  mathematical  model  is  given  in  Fig.  6. 

It  comprises  an  external  supersonic  diffuser,  a  short  cowl  and  an  ex¬ 
ternal  or  plug-type  nozzle,  all  two-dimensional  (or  possibly  half- 
round).  The  approaching  air  is  compressed  by  the  forward  isentropic 
ramp.  The  rim  of  the  cowl  is  placed  at  the  point  where  the  inlet  is 
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Fig.  6  MODIFIED  TWO-DIMENSIONAL  ENGINE  WITH  HIGH-SUBSOKIC 
HEAT  ADDITION,  "THERMAL  THROAT,"  AND  EXTERNAL 
SUPERSONIC  EXPANSION 
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designed  for  shock-on  rim  operation  for  a  given  flight  Mach  number. 
The  functions  of  the  cowl  are  to  capture  the  appropriate  quantity  of 
compressed  air.  to  separate  the  shock  wave  and  the  flame  zone  so 
that  there  is  sufficient  time  for  mixing  and  the  ignition  lag, and  to 
turn  the  burned  gas  to  the  desired  angle  0  over  the  rearward  surface 
so  that  the  supersonic  expansion  will  take  the  most  advantageous 
pattern  with  respect  to  the  required  thrust  and  lift  production. 

In  order  to  stabilize  the  flame  at  the  cowl  exit,  the  flow  under 
the  cowl  is  expanded  or  compressed  slightly,  depending  upon  the 
rate  of  heat  addition  (fuel-air  ratio),  which  in  turn  must  be  just 
sufficient  to  accelerate  the  gases  to  sonic  velocity  as  they  pass 
through  the  flame.  It  is  assumed  that  the  required  expansion  or  com¬ 
pression  is  accomplished  between  Stations  2  and  3,  and  that  the  fuel 
is  added  and  mixing  is  accomplished  in  a  constant  area  duct  between 
Stations  3  and  4.  Combustion  is  completed  at  Station  4  where  =  1. 

Preliminary  Assessment  of  Diffuser  Variables 
(Perfect  Gas  withy  =  1.4) 

Relative  effects  n  thrust  of  the  degree  of  compression  at  Sta¬ 
tion  1  and  of  subsequent  expansion  or  compression  to  various  M  's 

J 

under  the  cowl  (Table  I)  are  estimated  with  the  following  simplifying 
assumptions. 

1.  All  points  represent  critical  operation  with  shock  on  rim 

for  a  given  M^,  ~  ~  constant»  anc*  Mg  is  varied  by 

varying  Ay 

2.  Fuel  is  injected  in  a  constant  area  section  between  Sta¬ 
tions  3  and  4  (A„  =  A  )  with  negligible  momentum  in  the 
flow  direction,  so  that  the  momentum  balance  is 
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^  =  p4A4  ^  +  ^4M4  )  ~  ^3  P3A3  ^  +  y3M3  ‘ 

which  reduces  to 

p4  <  i  +  y4>  *  P3  U  +  y3M32)  (2) 

3.  All  gases  are  ideal  gases  with  a  constant  specific  heat 

ratio  (yQ  =  yg  =  y 4  =  =  1*4)‘ 

4.  The  burned  gas  undergoes  insentropic  expansion  (P  = 

*5 

P  )  from  Station  4  (P  determined  from  Eq.  (2); 

*4  4 

=  1)  to  Station  5,  where  it  is  fully  expanded  to 

ambient  pressure  (Pg  =  PQ;  M,.  obtained  from  resulting 

P,/ P.  for  y»  1.4). 
b  *5 

5.  The  cowl-body  geometry  is  adjusted  so  that  in  effect  a 
one-sided  isentropic  nozzle  is  provided;  that  is,  the  ef¬ 
fective  slip  plane  between  burned  gas  and  adjacent  air 
stream  is  parallel  to  the  flight  axis,  so  that  the  burned 
gas  flow  vector  at  Station  5  is  parallel  to  the  flight  axis, 
and  A,  =  A  .  (A  /A  is  calculated  from  M  . ) 

Or  D  4  D 

(Assumptions  3  and  5  are  admittedly  naive  but  are  con¬ 
sidered  adequate  to  the  present  purpose;  they  are  cor¬ 
rected  in  later  calculations. ) 

Thus  for  selected  values  of  P  ,  Mn,  M  ,  and  M_,  the  corre- 

U  U  1  0 

sponding  values  of  M,.  and  A  are  determined,  and  the  approximate 

D  0 

gross  thrust  is  calculated  from 

Fg  =  ^5  "  "  P0  (A5  '  A0) 
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which  for  convenience  is  nondimensionalized  and  re-written  in  the  form 


Wo 


■  <VAo> 


(1  +yM5  )  -  (1  +yMa  )  -  ((As/A0)  -  1)  (41 


a  thrust  coefficient  can  then  be  calculated  from 


V-F 


poAo 


Vo 

qOA5 


poAo 


■>'oM0  (A5/Ao’ 


(5) 


The  results  are  also  related  to  air  specific  impulse,  S  ,  a  parame- 

a 

ter  familiar  to  jet  engineers,  which  is  defined  as  the  critical  (M  =  1) 
stream  thrust  per  pound  of  air  per  second: 


S 

a 


P3A3(1  +  y3M32) 

g  poVo^Vo 


(Si- 


Table  I  gives  the  results  of  these  over-simplified  calculations 
for  a  range  of  M  's  for  both  M  =  2.  75,  which  represents  a  good  isen- 

w  A 

tropic  ramp  diffuser  (kinetic  energy  efficiency  of  0.  947  for  y  =  1.4),  and 
Mj  =  4,  which  lepresents  a  simple  10°  wedge.  The  following  general 
points  are  shown  by  this  table: 

1.  For  a  given  of  S  ,  the  thrust  coefficient  is  increased 

o  3 

by  60  percent  to  100  percent  when  is  decreased  from 
4  to  2.  75  by  better  diffusion. 

2.  The  above  statement  applies  for  the  range  of  M's  of 

O 

greatest  interest,  which  is  shown  to  be  the  range  of  0.  3 
to  0.5,  because  these  M ,'s  correspond  to  S  's  of  140 
to  200,  which  in  turn  correspond  to  the  values  obtained 
with  common  fuels.  For  example,  the  complete  com¬ 
bustion  of  kerosene  at  an  equivalence  ratio  of  0.5,  with 
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air  at  an  inlet  total  temperature  of  2360°F  provides 

an  S  of  160. 
a 

3.  The  above  ranges  in  M„  and  S  provide  a  2-1/2  - 

o  cl 

fold  variation  in  C,p  for  the  good  diffuser  (M^  =  2.  75). 

This  variation  is  obtained  by  a  fifty  percent  change 
in  the  nozzle  areas  A  and  A  .  These  numbers  point 

TT  D 

out  the  importance  of  variable  geometry  at  this  flight 
speed,  but  they  also  suggest  that  (as  is  usually  the  case 
for  ramjets)  higher  C  '  s  might  be  obtained  with  smaller 
A  _'s,  that  is,  without  full  expansion. 

D 

Effects  of  Operational  Variables  for  Fully 
Expanded  Real  Gas  Cases 

Effects  of  fuel-air  equivalence  ratio,  ambient  pressure,  flight 
Mach  number,  and  frozen  versus  equilibrium  exhaust  flow  were  de¬ 
termined  for  the  kerosene-air  system  (Ref.  22).  The  air  temperatures, 
pressures,  and  Mach  numbers  through  Station  3  were  calculated  as 
previously  except  that  caloric  imperfections  were  taken  into  account 
(Ref.  23).  Conditions  for  the  bur  ned  gas  at  Station  4  were  determined 
by  mass  balance,  Eq.  (7),  momentum  balance,  Eq.  (2),  and  energy 
balance,  Eq.  (8)  below,  between  Stations  3  anu  4,  allowing  for  the  fuel 
addition  between  these  two  stations  and  employing  the  assumption  that 

=  ag  ,  the  sonic  velocity  under  conditions  of  chemical  equilibrium 
4 

at  Station  4.  The  burned  gas  was  again  assumed  to  expand  isentropi- 
cally  to  P  =  P  ,  and  conditions  at  Station  5  were  calculated  for  both 

D  U 

frozen  and  dynamic  equilibrium  by  entropy  and  energy  balances  as 
follows. 
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Mass  balance: 


Energy  balance 


V3  (1  +  f)  *  >4  V4 


h.  +  fh, 

t  Y 

_A  T 

1  +  f 


V 


\  h4  +  2gJ 


V. 


ht5  =  h5  +  2gJ 


(7) 


(8) 


where  h  and  h  are  the  stagnation  enthalpies  of  fuel  and  air  respec- 
A  V 

tively,  and  other  enthalpies  are  per  unit  mass  of  burned  gas. 


Equations  (2),  (3),  and  (8)  are  satisfied  by  iteration.  One  route 

is  to  estimate  a  P4  (for  the  present  model  »  0.53  P^)  then  estimate 

the  T4  and  =  a  which  will  satisfy  Eq.  (8),  referring  to  tables  of 
4 

thermodynamic  properties  of  combustion  gases.  Now  these  values  of 
and  are  used  to  estimate  through  0(M^),  a  tabulated  function 
of  Mach  number  (see  Eq.  (9)  and  Ref.  24).  This  is  done  by  employing 
Eqs.  (2)  and  (7)  together  with  fhe  perfect  gas  law  and  the  relationships 
between  Mach  number  and  velocity  and  between  total  and  static  tempera¬ 
ture  (ignoring,  for  the  moment,  caloric  imperfections  which  could 
cause  0.2  percent  error  in  temperature  ratio)  to  show  that 


0(.)  a 

0(M3) 


(1  +  yM2)/(2(y  +  1)M2(1  +  M2)]1  ‘ 2 


R.T. 
,  4 

(1  *  f)-T7 


(1  T  v,  /R.T . 

4  4  4 


v4  nyn 

V2T7- 


r,  Vt. 


(9) 
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from  which  it  is  seen  that  0{M^)  is  primarily  dependent  on  T^,  V^,  and 

Tt  =  T  ,  because  appropriate  values  of  y„  and  will  fall  within 
t0  t  '  r3  4 

«j  u 

narrow  ranges.  Knowing  0(M,,)  from  Eq.  (9a)  one  could  then  solve  for 
M„  from  Eq.  (5)  or  by  interpolation  for  7/5  <  y  <  9/7  in  the  tables  of 
Ref.  24. 

The  pressure  P  and  area  A  =  A  may  obtained  from  normal 

O  v  * 

shock  tables,  corrected  for  caloric  imperfections;  is  then  recalcu¬ 
lated  from  the  momentum  balance  and  compared  with  the  assumed  value. 
The  process  is  iterated  until  the  difference  between  estimated  and  re¬ 
calculated  P^  is  too  small  to  noticeably  affect  and  the  corresponding 

V 

Exit  conditions  at  station  5  were  then  obtained  by  energy  (Eq.  8) 
and  entropy  balances.  For  dynamic  equilibrium,  in  which  the  chemical 
composition  is  allowed  to  shift  with  T  and  P,  the  properties  of  the 
burned  gas  may  be  read  directly  from  the  tables  of  Ref.  22  by  inter¬ 
polating  to  s.  =  s  at  the  exit  pressure  P  =  P  .  The  energy  balance 
0  4  0  u 

then  gives  V  .and  the  mass  balance  between  stations  4  and  5  gives  A_. 

0  0 

The  stream  thrust  a>.  station  5  is  then  calculated  from 

^5  =  fVV1+V5/R51V  <10) 


and  gross  thrust  is  calculated  from  Eq.  (3). 

For  frozen  equilibrium,  where  the  burned  gas  composition  (mol 
fractions,  mean  molecular  weight)  at  station  4  is  assumed  to  hold  for 
all  points  further  downstream,  the  equation  for  the  isentropic  process 
is 


2.X.  (S°/R).  -  ln{P  /P_ )  =  ^T.  X-  (S°/R). 

i  i4  Vj  4  0  i4  i5 
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where  the  mol  fraction  and  entropy  function  for  each  pure  product 

species  is  obtained  from  the  tables  at  T.,  P..  The  T.  which  satisfies 

4  4  5 

this  equation  at  P,.  =  P^  is  found  by  trial  and  error  from  interpolation 
in  the  tables.  The  specific  enthalpies  of  the  various  product  species 
at  this  temperature  are  then  used  to  calculate  the  burned  gas  specific 
enthalpy,  which  is  next  used  in  the  energy  equation  to  calculate  V  • 

R5  =  R4* 

Results  of  these  calculations  are  presented  in  '  ible  II.  Points 
to  be  noted  are 

1.  Thrust  coefficients  are  encouragingly  high  and  far  out¬ 
weigh  the  major  engine  drag  force,  the  cowl  lip  drag, 
which  is  crudely  estimated  at  the  bottom  of  the  table. 

2.  As  one  might  expect,  fuel  economy  as  represented  by 
fuel  specific  impulse  [lb  gross  thrust /(lb-fuel /sec)], 
decreases  slightly  as  fuel-air  equivalence  ratio  is  in¬ 
creased  from  0.25  to  0.70. 

3.  There  is  little  difference  between  frozen  and  dynamic 
equilibrium  results  at  Mach  5.  (Only  the  E.R.  =  0.  7  re¬ 
sults  are  shown  in  the  table,  but  there  was  even  less  dif¬ 
ference  at  lower  equivalence  ratios  as  one  might  expect. ) 

4.  Effect  of  altitude  in  the  isothermal  range  will  be  negligi¬ 
ble. 

5.  Thrust  coefficient  and  fuel  specific  impulse  fall  off 
rapidly  as  flight  Mach  number  approaches  8.  Since  drag 
would  decrease  only  very  slightly  as  speed  .increased,  a 
zero  net  chrust  condition  would  soon  be  reached  at  E.R.  = 
0.  5.  Note,  however,  that  the  overall  efficiency  of  the 
engine,  as  represented  by  the  product  of  velocity  and  fuel 
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Table  II 


Fully  Kxpanded  Operation  of  Isentropic  KRJ  with  Cov\l  effects  >f  Flijht  Speed, 
Pressure,  Fuel-Air  Ratio,  and  Chemical  Kquilibrij  : 

Real  Gas,  Kerosene-Air 


'•  UK 

KRJ 

RJ 

*  *  *  r  Mroam 

*0 

5 

5 

6 

8 

5 

V  :,im 

0.  1 

0.  10 

0.01 

0.01 

0.01 

,  °R 

2220 

30 

2950 

4740 

2220 

.  1  shock 

S 

2.  75 

2.75 

3.25 

4.  10 

2.75 

'  '2 

0.  1293 

0.  1293 

0.  0959 

0.  0499 

0.  1293 

0.483 

0.483 

0.447 

0.  325 

0.483 

1  .«<  1  I  quiv.  Ratio 

0.25 

0.50 

0.  70 

1.  0 

0.5  b 

0.50 

.50 

0.50 

i'i  .>”,u  Diff. 

V,J 

0.488 

0.393 

0. 35! 

0.346 

0.392 

0.445 

C.  561 

0.  044 

P  ,  atm 

1.919 

2.  021 

2.060 

2.066 

20.  02 

3.920 

11.  10 

2.236 

2135 

2167 

2175 

2180 

2.  : 

2860 

4525 

2219 

'•  -  '  A4 

0.  1283 

0.  1519 

0.  1670 

0.  1694 

0.  153 

0.096 

0.  032 

1.250 

.1-0  it,  M  ,  =  1* 

4 

V  fps 

2488 

2761 

2907 

3012 

2733 

2926 

3240 

2765 

P,,  atm 

1.  100 

1.082 

1.087 

1.112 

10.  79 

2.218 

7.  13 

1.234 

.,  ’R 

2786 

3542 

4058 

4469 

3555 

4059 

5107 

355 1 

.  'em.  Comp.  ’■« 

D 

D 

F 

D 

D 

D 

D 

D 

D 

S'  fp< 

5760 

6580 

7033 

7120 

7765 

6570 

7370 

884  0 

6410 

1  . ,  R 

862 

1201 

1428 

1485 

1920 

1204 

1200 

1291 

1425 

1.877 

2.320 

2.  640 

2.698 

2.040 

2.337 

2.490 

2.  900 

i.  250 

1  _  'PA 
.  0  0 

43.  98 

51.  22 

55.84 

56.35 

58.54 

51.  17 

68.  07 

107.9 

50.  24 

Pi  rformance 

CT  "  Fg/q0AF 

o.  3 16 

0.342 

0.394 

0.395 

0.385 

0.338 

0.212 

9.  1 19 

0.639 

• 

*f  Fg/Wf 

1810 

1780 

1660 

1700 

1310 

1770 

1C20 

1240 

1790 

Drag,  C  *** 

10°lip 

0.  053 

0.  065 

0.072 

0.  072 

0.  077 

0.  066 

0.064 

0.062 

0. 023 

.stagnation  specific  enthalpy  of  air  =  565  Btu/lb  at  2220°R,  fuel  assumed  to  enter  at  i77°F  with 
•189  Btu/lb. 

!  -  Composition  frozen  beyond  station  4.  D  =  composition  shifting  with  T  and  P  in  dynamic  equil. 
Wave  drag  of  10°  cowl  Up  on  projected  area  Aj,-A^. 
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specific  impulse  increases  with  Mach  number  (it  may 
pass  through  a  maximum  between  6  and  8). 

6.  Comparison  of  the  last  column  for  a  ducted  ramjet  with 

the  second  column  for  the  ERJ,  both  with  the  same  assump¬ 
tions  on  diffuser  total  pressure  recovery  (39.4  percent) 
and  combustion  efficiency  (100  percent),  shows  that  the 
ERJ,  if  it  could  in  fact  obtain  full  isentropic  expansion, 
would  produce  comparable  1^.,  but  would  have  a  much  lower 
and  higher  because  of  its  large  frontal  area.  The 
arbitrarily  chosen  area  ratios  for  the  ramjet,  which  are 
typical  of  current  practice,  result  in  under  expans  ion.  Sub¬ 
sequent  calculations  will  show  that  the  ERJ  also  benefits 
from  a  reduction  in  drag  at  the  penalty  of  underexpansion. 

Effects  of  Cowl- Body  Geometry,  Varying  Expansion 
Real  Gas  Cases 


A  few  cases  of  underex  pans  ion  were  solved  to  obtain  a  feeling  for 
the  balance  between  the  effects  of  cowl  drag  and  degree  of  expansion  on 
the  gross  thrust  and  lift  of  the  system.  The  method  of  characteristics 
was  used  to  obtain  the  pressure  distribution  on  the  surface  aft  of  the 
cowl. 

The  same  assumptions  regarding  the  approach  flow  and  plane 
flame  zone  are  made,  so  that  the  conditions  at  station  4,  immediately 
downstream  of  the  flame  front,  are  the  same  as  those  given  in  Table  II 
for  Mq  =  5,  Pq  -  0.  01  atmosphere,  and  the  appropriate  equivalence 
ratio.  For  simplicity,  it  is  assumed  that  the  after  surface  is  an  isen¬ 
tropic  surface  and  that  Prandtl-Meyer  expansion  occurs  from  the  cowl 
lip.  For  a  given  A  / An,  the  ratio  of  cowl  escape  area  to  capture  area. 
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and  a  given  flow  angle  0^  on  the  downstream  side  of  the  flame,  the 
parameters  required  to  calculate  the  forces  at  the  last  Mach  line  to  hit 
the  missile  are  determined  as  follows.  The  pressure  and  temperature 
at  any  Mach  line  are  determined  as  functions  of  V  by  energy  and 
entropy  balances,  assuming  dynamic  equilibrium  in  the  burned  gas 
(see  Eq.  (8)  and  previous  discussion  on  dynamic  equilibrium).  From 
the  tables  of  properties  of  combustion  gases,  the  acoustic  velocity  a, 
and  hence  the  Mach  angle,  $  =  sin  *  a/V,  may  be  determined  as  a 
function  of  V.  The  local  flow  angle  0  is  then  determined  from 


V 


dV 


0  =  0  -  J  cot[3~  . 


(12) 


I 

1 

I 

I 


The  area  A,.,  represented  by  the  length  of  the  last  Mach  line,  is  deter¬ 
mined  by  mass  balance 


A5P5  V5  s*n^5  "  P4Y4A4 


(13) 


A  is  related  to  Aj.  by 
6  0 


Ae  ■  A5  sin  «5  +  J5) . 


The  stream  thrust  across  the  last  Mach  line  is 


(14) 


3*5  -  P5Ae+(PSASV5StoVR5T5,<V5cos95>  '  <I5) 


These  relationships  may  be  seen  from  Fig.  7. 

Using  this,^,  thrust  coefficient  and  fuel  specific  impulse  are 
calculated  as  before.  Lift  forces  are  calculated  for  the  same  force 
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Fig.  7  MODEL  FOR  CALCULATION  OF  EXIT  FORCES  AT  LAST 

MACH  LINE  STRIKING  NONREFLECTING  EFFUSOR  SURFACE 

envelope  used  in  the  thrust  calculation.  Since  the  bounding  streamline 
between  stations  0  and  1  is  parallel  to  the  flight  axis  and  undisturbed, 
no  momentum  is  transported  across  it,  and  the  pressure  force  on  it  is 
equal  and  opposite  to  the  pressure  force  on  the  corresponding  portion 
of  the  flat  top  missile  surface.  Hence  the  left  forces  need  be  evaluated 
only  for  the  last  Mach  line,  the  cowl,  and  the  corresponding  areas  on 
the  top  surface. 

L5  =  (P5  -  P0>  A5  COS  <V  +  <P5  Vo  Si"  VRo  V  • 

ret 

(V.  sin  0_), 


(16 
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where  A  is  the  longitudinal  projected  area  of  the  cowl,  and 

CL  =  (L,  +  L  )/q0A  .  (13) 

net  net 

(Further  details  for  the  calculation  of  L  are  given  in  Appendix  I. ) 

c 

For  any  chosen  0  (flow  angle  leaving  the  cowl)  and  A  /A  (ratio 
4  0  u 

of  cowl  escape  ares  to  capture  area),  one  can  calculate  C^,  and  1^.  from 

the  gross  thrust  and  CT  from  Eq.  (18).  For  example,  one  can  choose 

.L/ 

a0,,  sav,  one  radian,  and  varv  A  /A,_  to  obtain  curves  of  the  above 
4  e  U 

parameters  versus  A^/A^.  As  A^/A^is  increased,  one  will  approach 
the  point  in  the  calculations  where  the  last  Mach  line  to  strike  the  mis¬ 
sile  will  have  turned  the  flow  parallel  to  the  flight  axis;  that  is,  the 
velocity  vector  Vg  will  be  parallel  to  the  flight  axis  and  0,.  will  be  zero. 
Further  isentropic  expansion  without  reflection  of  Mach  lines  would 
require  the  expansion  surface  to  turn  away  from  the  flight  axis  Egain, 
a  case  which  would  probably  not  be  desirable  in  the  overall  missile 
design.  (It  would  result  in  decreased  thrust  and  increased  drag  and 
would  perhaps  make  the  engine  more  difficult  to  integrate  with  an  air¬ 
frame;  it  would,  however,  provide  more  lift,  but  a  fiat  plate  extension 

•'•ould  also  accomplish  this  without  the  thrust  penalty, )  At  any  rate,  it 

2 

was  decided  to  terminate  these  calculations  at  0C  =  0. 

o 

Figure  8  illustrates  the  general  configurations  of  tne  sets  of 
curves  which  can  be  developed  in  this  manner  and  the  way  that  they  are 

2  ‘  ~~ 
One  might  ask  whether  it  is  valid  to  carry  these  calculations  for  a 

Prandtl-Meyer  expansion  even  to  the  point  0^  =  0.  Examples  in 

Appendix  II  show  that,  at  least  in  the  region  of  greatest  interest,  this 

procedure  is  valid. 
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Fig.  8  EXAMPLE  OF  PERFORMANCE  MAPS  FOR  EERJ  PRODUCED 
BY  VARYING  Aq/aq  AT  VARIOUS  S4's 

Curves  are  terminated  at  «5  -  0,  which  is  consid¬ 
ered  the  limit  of  practical  interest.  Numerical 
values  are  not  comparable  to  subsequent  figures 
because  they  are  incorrectly  based  on  a. 
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terminated  at  the  "locus  of  9r  =  0. "  This  figure  is  presented  here  for 

o 

two  reasons:  to  illustrate  the  shapes  of  the  curves,  and  to  point  out  the 
considerable  erroi  s  which  can  be  introduced  into  such  calculations  by  the 
use  of  improper  acoustic  velocities.  In  calculating  the  curves  for  Fig.  8, 
it  was  assumed  that  the  burned  gases  were  in  dynamic  equilibrium;  how¬ 
ever,  it  was  erronec  usly  assumed  that  since  a  characteristic  (Mach  wave) 
represents  an  infinitely  t*  in  mathematical  surface,  there  would  be  no 
time  for  chemical  change  in  crossing  it  and  the  frozen  sonic  velocity 
should  be  used  even  in  dynamic  equilibrium  calculations.  This  is  rot 
correct,  because  in  the  basic  potential  flow  equation  from  which  differ¬ 
ential  equations  are  obtained  for  the  characteristics,  the  local  value  of 
a  is  used,  and  this  value  is  determined  by  the  thermodynamic  nature  of 
the  gas,  that  is,  a  =  dp/d p.  Therefore  in  dynamic  equilibrium  calcula¬ 
tions  one  should  use  1  o  equilibrium  acoustic  velocity  a ^  rather  than 
the  "frozen"  acoustic  velocity  af. 

The  numerical  values  of  Fig.  8  should  not  be  used  by  the  reader 
for  any  further  purpose,  because  the  C^,  and  1^  values  are  erroneously 
high,  whereas  the  lift  coefficients  are  low. 

Rather  than  recalculate  all  of  the  curves  generated  for  figures 
such  as  Fig.  8,  it  was  decided  that  an  investigation  in  greater  detail 
be  made  of  the  more  inte  esting  but  restricted  case  which  represents 
the  maxima  in  var;ous  parameters  for  isentropic  expansion  without 
wall  reflections,  namely,  the  case  of  0,.  =  0,  using  the  correct  equili¬ 
brium  sound  velocity.  This  was  done  by  setting  0  -  0,.  =  0  in  Eq.  (12)  so 
that  for  each  0^  there  is  a  corresponding  value  of  the  local  velocity 
which  will  have  been  reached  when  the  flow  has  become  parallel  to  the 
flight  axis  (or  flat  top  of  the  engine). 
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The  basic  parameter  being  investigated  in  this  case  is  the  degree 
of  under  expans  ion,  that  is,  the  ratio  PK /P  at  the  point  where  0  =  0. 

Jn  other  words,  the  cowl-body  geometry  is  varied  to  produce  various 
degrees  of  under  expans  ion,  but  in  all  cases  it  is  a  turned  that  the  geom¬ 
etry  is  designed  to  give  the  maximum  (0  =  0)  value  of  the  performance 

D 

parameter  corresponding  to  that  P.  /P„.  The  iirst  figure  in  this  series, 

5  O 

Fig.  9,  therefore  pi  esents  the  relat’onshiDs  between  Pc  / P^  and  the 

b  (j 

geometry  variables,  0^  and  A^/A^,  required  to  produce  =  0  at  vari¬ 
ous  kerosene-air  equivalence  ratios  for  a  flight  speed  of  Mach  5  and  P^ 
0.01  atmosphere  and  T^  -  395°R.  The  results  are  thereafter  presented 

as  functions  of  P  /P  and  are  discussed  in  the  following  sec  .ion. 

0  \) 

Discussion  of  Results  for  Varying  Expansion.  Real  Gas 

Figure  10  presents  the  thrust  coefficients  for  various  cases  as 
functions  of  the  degree  of  underexpansion,  P^/P^,  with  geometry  ad¬ 
justed  to  give  0  =  0,  as  discussed  in  the  previous  section.  All  of  the 

D 

calculations  assume  100  percent  combustion  efficiency..  The  uppermost 
curve  in  Fig.  10  also  assumes  a  nozzle  velocity  coefficient  of  unity  (no 
losses  in  expansion).  All  other'  curves  in  this  and  subsequent  figures 
are  based  on  a  nozzle  velocity  coefficient  C,.  of  0.  97.  Here  it  is 
assumed  that  any  deficiency  in  mass  flow  which  would  have  resulted 
with  a  theoretical  nozzle  (boundary  layer  effects)  has  been  compensated 
in  the  nozz"  e  design,  so  that  the  0.  97  coefficient  applies  only  as  a 
first  power  efiect  on  velocity  V_. 

Effect  of  Nozzle  Velocity  Coefficient  on  C _ 

The  difference  between  the  first  and  second  uppermost  curves 
of  Fig.  10  indicates  the  effect  of  a  3  percent  loss  in  exit  velocity  when 
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Fig.  1C  THRUST  COEFFICIENTS  VERSUS  EXPANSION  RATIO, 
p5/p0,  FOR  MODIFIED  ENGINE  WITH  eg  -  0, 

AT  VARIOUS  EQUIVALENCE  RATIOS 

Upper  curve  for  E.  R.  -  1-0  based  on  nozzle 
velocity  coefficient  C  -  1.00;  solid  curves 
based  on  C  «  0.97  (see  text);  dashed  curves 
represent  net  thrust  coefficient  (C  -  0-97) 
after  subtracting  estimated  drag.  Sther 
assumptions  arc  the  same  as  in  Fig.  9- 
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burning  kerosene  at  an  E.R.  of  1.  0.  The  thrust  loss  is  approximately 
9  percent.  Each  point  lost  in  C^.  therefore  results  in  a  3  percent  loss  in 
thrust  at  E.R.  =  1.  0,  and  could  result  in  as  much  as  6  percent  loss  at 
this  equivalence  ratio  if  the  geometry  were  not  corrected  for  boundary 
layer  effects  or  were  operated  off-design  (for  example,  by  simply 
adjusting  the  cowl  over  a  body  that  was  designed  for  different  operating 
conditions)..  At  lower  equivalence  ratios  the  effect  of  on  C^,  is  even 
greater,  because  the  difference  between  exit  stream  thrust  and  the  fixed 
entrance  stream  thrust  decreases  with  equivalence  ratio,  so  that  a  given 
loss  in  exit  stream  thrust  is  magnified  in  its  effect  on  gross  thrust  and 
thrust  coefficient.  For  example,  at  E.R.  =  0.25  each  point  lost  in 
would  cause  a  6  percent  loss  in  maximum  C^,  (with  a  mass-flow  com¬ 
pensated  geometry). 

Effe  .  of  Under  expans  ion  on 

All  of  the  curves  of  Fig  10  illustrate  that  better  thrust  coeffi¬ 
cients  (based  on  required  engine  frontal  area)  are  indeed  obtained  when 
t.ie  burned  gas  is  underexpanded  as  it  leaves  the  isentropic:  expansion 
surface.  For  the  richer  mixtures,  E.R.  =  0.  7  to  1.0,  the  thrust  coeffi¬ 
cients  peak  rather  sharply  at  an  exit-to-ambient  pressure  ratio  of 
seven  to  eight.  For  lower  PK  /P~.,  the  frontal  area  A  must  be  in- 
creased  in  order  to  provide  the  A^lA^  corresponding  to  the  flow  turn¬ 
ing  angle  0  which  is  necessary  to  reach  that  P  /P  •  since  A  in- 
creases  more  rapidly  chan  the  gross  thrust  (with  the  present  assumptions 
on  geometry.  Appendix  I)  C.^,  based  on  frontal  area  decreases. 

At  lower  equivalence  ratios,  the  peak  occurs  at  lower  P  /P 

D  U 

but  even  at  E.R.  =  0.25,  which  is  probably  near  the  practical  lower 
limit  with  regard  to  gross  thrust  fuel  economy  (Fig.  11),  a  pressure 

ratio  Pc/P^  =  3  should  be  maintained. 

o  (J 


-  37  - 


tftiinre'HiftBfimu 


a 


-  hf  • 


TMt  JOHN*  KONKIN*  UNIVCftSITY 

APPLILD  PHYSICS  LABORATORY 

*U.VKH  SM)M  maktvano 

Effect  of  Equivalence  Ratio  on  C^. 

The  solid  curves  of  Fig.  10  compare  C^'s  with  =  0.  97  ,  for 
equivalence  ratios  of  1.0,  0.7,  0.5,  and  0.25.  Maximum  CT  is  roughly 
proportional  (and  early  equal  to)  equivalence  ratios  in  the  range  of  0.25 
to  0.  7,  but  increases  more  slowly  from  0.  7  to  1.  0  and  probably  falls 
off  beyond  E.R.  =  1.0.  The  latter  effect  is  more  pronounced  on  either 
side  of  the  optimum  expansion  ratio,  that  is  the  E.R.  =0.7  and  E.R.  = 
1.0  curves  converge  for  either  high  or  low  expansion. 


i 

1 


t 


Effect  of  Engine  Drag 

The  difference  between  the  solid  curve  and  the  dashed  curve 
for  each  equivalence  ratio  in  Fig.  10  represents  the  estimated  engine 
drag  coefficient.  (See  Appendix  I  for  method  of  estimation. )  Pressure 
ratios  of  5  to  50  represent  short  engines  with  relatively  low  friction 
drag  and  cowl  pressure  drag,  so  that  the  estimated  engine  drag  coeffi¬ 
cients  fall  in  the  range  0.  05  to  0.  07.  For  greater  or  lesser  expansion, 

P  / P  <  5  or  <  50,  the  engine  must  be  lengthened,  and  the  estimated 
5  0 

drag  coefficients  increase  to  0.  09  to  0.  11.  Moreover  the  latter  values 

probably  should  be  even  higher,  because  the  geometry  assumptions 

lead  to  a  long  subsonic  duct  for  either  very  large  or  very  small  A  / A^ 

“60 

(Fig.  9),  and  no  allowance  was  made  in  these  drag  estimates  for  in¬ 
ternal  friction  losses. 


Fuel  Specific  Impulse 


It  is  not  surprising  that  the  maximum  (defined  as  gross  thrust 

per  pound  of  fuel  burned  per  second,  I.  =  F  / w.)  is  obtained  with  com- 

f  g  f 

plete  expansion  in  all  cases.  Figure  11  indicates  that  the  best  fuel 
ec  >nomy  for  this  engine,  burning  kerosene,  will  occur  with  an  equiva¬ 
lence  ratio  between  0.5  and  0.  7.  At  maximum  conditions,  with 
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E.R.  ~  0.  6  and  P_ /P^~  7,  a  maximum  practical  value  of  1300  and  1400 
5  O 

appears  to  be  appropriate.  The  following  section  will  show  that  maximum 

range  might  be  obtained  with  slightly  more  expansion,  P  /P  ~  3,  where 

5  O 

an  linear  1500  would  be  appropriate. 

Optimum  Expansion  for  Lift 

Lift  coefficients  are  shown  in  Fig.  12,  referenced  to  engine 
planform  area,  CT  ,  read  from  solid  curves  and  left-hand  ordinate,  and 
to  engine  frontal  area,  CT  ..  ,  read  from  dashed  curves  and  right-hand 

ordinate.  These  lift  coefficients  peak  at  higher  P^/P^  than  do  the 
thrust  coefficients,  because  the  ratio  of  lift  to  thrust  increases  as  the 
initial  flow  angle  leaving  the  cowl,  0^,  is  decreased.  In  further  explana¬ 
tion,  the  lift  and  thrust  forces  on  the  rear  of  the  engine  could  be  obtained 
in  a  different  way  by  adding  their  respective  components  owing  to  the 
momentum  of  the  burned  gas  leaving  the  cowl  at  station  4  and  owing  to 
the  integrated  pressure  force  over  the  expansion  surface  between  station 
4  and  the  trailing  edge  (refer  to  Fig.  7).  Since  the  velocity  vector 
has  an  upward  component,  the  stream  thrust  which  produces  an  equal 
and  opposite  force  on  the  engine  at  this  point  has  a  downward  component, 
hencetj^  contributes  negatively  to  lift  according  to  sin  0^  and  positively 
to  thrust  in  proportion  to  cos  0^.  Moreover,  as  0^  is  decreased,  the 
lift  component  on  the  effusor  surface  owing  to  expansion  of  the  burned  gas 
will  increase  relative  to  the  thrust  component. 

The  lift  coefficients  must  eventually  pass  through  a  maximum, 
with  the  present  assumptions  on  cowl  geometry  (Appendix  I),  because 
both  planform  and  frontal  area  increased  area  increase  rapidly  for 
Pj. /Pq  >  10,  and  in  this  region  the  cowl  contribution  to  lift  is  an  in- 
c  easing  negative  quantitv.  On  the  other  hand,  for  P^/P^<  5,  the 
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Fig.  12  LIFT  COEFFICIENTS  VERSUS  Pg/Pg  FOR  VARIOUS 
EQUIVALENCE  RATIO.6 

CL  based  on  engine  planforn  area;  based  on 

engine  frontal  area.  Conpare  CL(AF)  to  CT  N£T  for 
relative  lift  and  net  thru3t  forces. 
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cowl  presents  an  increasing  positive  contribution  to  the  lift  force,  a 
case  analogous  to  a  conventional  engine  plus  wing. 

Relative  lift  forces,  which  include  the  ccwl  contribution,  should 
be  compared  to  net  thrust  forces  through  C^/A  ^  and  C^,  N£rp- 

Effect  of  Equivalence  Ratio  onC^ 

For  an  E.  R.  of  s  0.  7,  the  effect  of  equivalence  ratio  on  CT  is 

i-t 

smaller  tnan  on  C_,  because  as  equivalence  ratio  is  de- 

I ,  JNii/ 1 

creased,  the  velocity  decreases  (see  Table  II),  and  hence  the 
adve:  sej^  effect  on  lift  discussed  in  preceding  paragraphs  becomes 
smaller. 

Lift  /Drag  Ratio  for  the.  Engine 

If  this  engine  could  be  stabilised  more  easily,  particularly 
against  the  pitching  moment  resulting  from  the  lilt  under  the  effusor 
surface  (unless  the  c.g.  is  far  to  the  rear),  then  it  would  provide  a 
vehicle  of  very  high  L ID  indeed.  Figure  13  shows  that  a  value  of  15  is 
estimated  for  stoichiometric  burning  with  Pj. /P^  =  8.  An  L/D  of  9  is 
estimated  even  for  the  lowest  equivalence  ratio,  for  which  the  corre¬ 
sponding  net  thrust  coefficient  would  be  only  (X  1.  However,  these  high 
L/D's  could  only  apply  for  climb  or  acceleration;  for  constant  velocity 
cruise  the  net  thrust  would  have  to  be  opposed  by  an  added  external 
drag. 

Optimization  for  Range 

The  significance  of  the  results  of  the  foregoing  calculations 
with  regard  to  a  long-range  cruise  vehicle  can  be  better  understood  by 
considering  the  Breguet  equation  for  maximum  range  at  constant 
velocity. 
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R  =  V(If  L/D)  ln(W0/W1) 


In  this  equation  thrust  equals  drag  and  lift  equals  weight.  For  maximum 
range,  one  wants  to  maximize  the  quantity  1^  L/D.  (If  changes  in  1^.  L/D 
noticeably  affect  the  logarithm  of  the  initial  to  final  weight  rttio  Wq/W^, 
then  one  must  maximize  (I^L/D)  In  (W^/W^),  but  the  latter  refinement 
is  beyond  the  scope  of  the  present  report.  Centrifugal  lift,  which 
amounts  to  about  3.  7  percent  of  vehicle  weight  at  Mach  5,  is  also 
neglected. ) 

As  a  first  approximation  to  a  vehicle  fulfilling  thr  requirement 
of  drag  equal  to  thrust,  it  is  assumed  that  the  most  economical  way  to 
increase  drag  is  to  add  wings.  An  L/D  of  6  for  added  wings  is  con¬ 
sidered  appropriate  for  Mach  5  cruise.  The  wing  drag  coefficient  based 
on  engine  frontal  area  is  assumed  to  be  equal  to  the  net  thrust  coefficient 
C  .  The  total  lift  coefficient  based  on  engine  frontal  area  is  then 

JL  *  JN  Hi  X 

taken  to  be  the  engine  lift  coefficient  ^  plus  six  times  CT 

F 

The  resulting  L/D  for  the  missile  is  this  total  lift  coefficient  divided 
by  the  gross  thrust  coefficient  (or  total  drag  coefficient)  C^. 

Figure  14  presents  the  total  lift-drag  ratios  estimated  in  this  way. 

An  overall  L/D  of  6.  0  to  6.  5  in  the  region  of  interest  for  E.  R.  =  0.  5  *o 

1.  0  is  indicated.  (The  E.R.  =0.25  curve  runs  off  the  plot  to  some 

maximum  L/D  near  12  at  P  /P  30;  values  for  higher  P_  / P  are  not 

0  U  «)  u 

of  interest  because  1^  is  negative. )  These  values  may  be  compared  wich 
the  L/D  range  of  4  to  5  often  mentioned  for  a  Mach  5  wing-body  vehicle 
with  a  conventional  ramjet  to  show  that  a  20  to  60  percent  gain  in  L/D 
might  be  expected  for  a  cruise  vehicle  with  external  burning. 
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Fig.  14  OVER-ALL  LIFT/DRAG  RATIO  FOR  CRUISE  VEHICLE  WITH 
ZERO  NET  THRUST,  BASED  ON  ADDITION  OF  WINGS  OF  L/D  -  6 
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Figure  15  shows  the  range  parameter,  1^  L/D  ,  for  the  engine  - 
plus-wing  vehicle  using  the  above  overall  L/D.  For  each  equivalence 
ratio,  the  maximum  L/D  occurs  at  P  /P  ~  3;  this  optimum  pres- 

x  D  U 

sure  ratio  becomes  more  sharply  defined  as  the  equivalence  ratio  is 
decreased.  The  highest  1^  L/D  calculated  was  8850  for  E.R.  =0.5  and 
Pj./Pq  =  3,  with  values  for  E.R.  =  0.7,  0.25,  and  1.  0  following  in 
decreasing  order.  Apparently  the  optimum  equivalence  ratio  would  be 
near  0.6  and  the  corresponding  L/D  would  be  9000  ±  100. 

On  the  basis  of  these  calculations,  an  engine-plus-wing  vehicle 
designed  for  cruise  at  Mach  5  at  100, 000  feet  with  an  equivalence  ratio 
of  0.5  and  P^/P^  =  3  would  have  the  following  characteristics: 


Minimum  ratio  of  engine  length  to  height,  i t  / h^ 
Corresponding  minimum  engine  body  volume,  vol/A 

r 

Percentage  of  lift  contributed  by  engine 

W  /A  at  start  of  cruise  =(CT  .  ,  .)  q_.A  /A_ 

O  F  2  '  L,  total  X)  plan  F 

For  50  I'd /ft  wing  loading,  wing  area  ^,/Ay 
Lift /drag  ratio  in  cruise 
Ij. L/D  in  cruise 


=  7.8 

3  ,  2 
=  3.  7  ft  /ft 

=  18  percent 

=  1000  lb  /ft2 

=  20  ft2 /ft2 

=  6.0 

=  8850  lb-sec /It 


Comparable  values  for  the  conventional  ramjet  of  Table  II  with 
Cy  =  0.  97  and  with  added  wings  of  L/D  =  6  would  be 

If  (with  Cv  =  0.  97)  1600  lb- sec /lb 

L/D  in  cruise  5.05 

L/D  in  cruise  8080  lb- sec /lb 

Thus  if  Wq/ were  comparable  for  the  two  engines,  the  range 
for  the  conventional  ramjet  vehicle  would  be  8  percent  smaller.  A 
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Fig.  15  OVER-ALL  RANGE  PARAMETER,  I-  L/D,  FROM  FIGS 

11  AND  14  1 

Included  for  comparison  at  E.  R.  -  0.5 
is  curve  for  effusor  rotated  downward 
through  65  degrees  (last  Mach  line 

parallel  to  flight  axis,  0 
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comparison  of  Wq/ W.  is  not  readily  made,  because  on  the  one  hand  the  ERJ 
requires  less  added  wing  area  per  square  foot  of  frontal  area,  less  cooling 
water,  and  less  nozzle  surface  area,  but  on  the  other  hand  two -dimensional 
engine  construction  would  be  heavier  than  axisymmetric  construction. 

One  can  see  that  if  the  engine  body  volume  of  the  above  ERJ 

vehicle  wer~  to  represent  the  fuel  storage  and  packaging  space  for  the 

missile,  the  minimum- length  engine  could  carry  very  little  fuel  per 

square  foot  of  frontal  area  and  would  have  a  range  of  the  order  of  500  to 

1200  nautical  miles.  The  above  numbers  would  therefore  te  appropriate 

only  to  small  short-range  missiles.  However,  if  the  body  under  the  cowl 

is  simply  lengthened  sufficiently  to  give  an  average  fueled  missile  den- 
3 

sity  of  80  lb /ft  at  the  beginning  of  cruise,  then  a  new  external  drag  coef¬ 
ficient  applies,  and  the  following  characteristics  are  estimated. 


Engine  length  of  height  ratio,  *t/hf 
Engine  body  volume,  vol/A 

Jr 

Lift /Drag  ratio  in  cruise 
If  L/D 

Cruise  fuel,  occupying  80  percent  of  body  volume 

at  50  lb/ft^,  w./A_ 
f  F 

Cruise  range 


=  14.2  ft /ft 
=  12.5  ft3 /ft2 
=  a.  73 

=  8420  lb-sec /lb 

=  500  lb /ft2 
=  4700  n.  miles 


As  noted  previously,  internal  friction  loss  under  the  cowl  is 

neglected.  A  rough  estimate  indicated  P  / P  ~  0.958,  or  AP  /q~  0.39, 

l3  l2  1 

ior  the  latter  example.  If  one  assumes  that  to  a  first  approximation  the 

exit  thrust  components  are  reduced  in  ratio  to  P  ,  then  1  is  reduced 

*3  f 

from  1470  to  1230  lb-sec /lb,  and  the  range  is  reduced  from  4700  n.  miles 
to  3900  n.  miles.  This  would  still  be  an  admirable  cruise  lange  for  a 
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Mach  5  missile  with,  say,  an  engine  1.  4  feet  thick  by  4  feet  wide  by 

2 

20  feet  long,  with  110  ft  of  added  wings  and  with  a  dry  weight  of  2800 

3 

pounds  and  a  fuel  capacity  of  2800  pounds.  Some  10  to  14  ft  of  body 
volume  would  remain  for  payload  and  controls;  this  payload  volume 
could  be  increased  by  carrying  part  of  the  fuel  in  the  wings.  Edge 
losses  for  the  engine  might  be  compensated  by  favorable  wing-body  in¬ 
terference  effects  in  an  integrated  design  (see  also  Section  V).  Use  of 
a  high  energy  fuel  of  25  percent  higher  heat  of  combustion  and  10  percent 
greater  density  would  increase  range  30  percent,  raising  the  3900  n.  miles 
figure  cited  above  to  5100  n.  miles. 


Effect  of  Vectoring  Nozzle 


For  any  given  equivalence  ratio  and  P_/P_,  the  lift  coefficient 

5  U 

can  be  increased  at  the  expense  of  thrust  coefficient  by  vectoring  the 
after  end  of  the  engine  downward.  A  few  calculations  were  made  to  in¬ 
vestigate  the  effect  of  vectoring  on  the  overall  1^  L/D.  It  was  assumed 
that  the  isentropic  external  diffuser  should  be  held  at  zero  angle  of 
attack  for  greatest  efficiency,  while  the  after  end  of  the  engine  was 
rotated  downward  to  put  Vj.  a.  a  positive  angle  of  attack.  However,  for 
the  cases  calculated,  the  loss  in  I  was  always  equal  to  or  greater  than 
the  gain  in  overall  L/D  (with  added  wings)  for  that  expansion  ratio  pro¬ 
viding  the  highest  oroduct  of  1^.  and  L/D.  In  other  words,  the  brief 
study  made  revealed  no  way  to  improve  on  the  maximum  overall  1^ 

L/D  which  was  calculated  for  the  straight  case  with  parallel  to  the 
flight  axis. 

Figure  15  includes  an  1^.  L/D  curve  calculated  for  an  engine  with 
the  after  end  rotated  downward  through  0^  degrees,  so  that  for  each 
point  on  the  curve  the  last  Macn  line  is  parallel  to  the  axis  and 
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makes  an  angle  5^  with  the  flight  axis.  The  top  surface  of  the  engine 

is  assumed  to  be  ’boattailed "  at  the  angle  flc.  (Refer  to  sketch  at  top  of 

o 

Fig.  16. )  Only  for  P^/Fq  near  unity  does  the  vectored- effusor  engine 

provide  reasonable  I.  L/D,  and  even  at  P,.  /P_  =  1,  the  value  produced 

i  a  o 

is  only  87  percent  of  the  maximum  1^  L/D  for  the  "straight  engine.'1 
Some  improvement  for  the  vectored  engine  could  be  obtained  by  stream¬ 
lining  the  top  surface  at  the  expense  of  body  packaging  volume,  but  it 
is  doubtful  that  it  would  produce  higher  L/D  than  the  straight  engine. 


Comparison  with  a  Conventional  Ramiet  with  Vectored  Nozzle 


On  the  basis  of  gross  thrust  versus  lift  force  due  to  internal 

flow  (that  is,  neglecting  lift  due  to  the  cowl),  the  conventional  ramjet 

and  the  external  expansion  ramjet,  both  with  vectored  nozzles,  give 

nearly  the  same  result  for  the  same  nozzle  velocity  coefficient  and  the 

same  exit-to-ambient  pressure  ratio  Pg/P^.  This  is  illustrated  in  Fig. 

16  in  which  F  /aA  is  plotted  both  against  L_  /  q„  A^,  where  L_ 

g  X)  O  5  O  u  5 

net  net 

is  "gage  corrected"  by  subtracting  the  product  of  Pq  and  thu  horizontal 

projection  of  for  each  engine,  and  against  /q^A^  where  the  "gage 

correction"  has  not  been  applied.  The  L_  curve  is  essentially  the 

5net 

same  for  both  engines;  the  small  advantage  for  the  conventional  engine 
owing  to  its  lower  total  pressure  loss  for  heat  addition  at  lower  sub¬ 
sonic  Mach  number  (Table  II)  appfa-ently  is  compensated  with  the  ERJ 
by  the  fact  that  the  gage  pressure  x  area  term  in  the  lift  force  is  greater. 
The  difference  between  L  and  Lj.  is  greater  for  the  ERJ  for  the  same 

D  D 

,,  ,  .  net 

reason,  that  is, 

Ag  cos  (8v  -  SgJP^gage  correction  for  ERJ 
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Fig.  16  VECTORED-NOZZLE  RAMJET  VERSUS  VECTORFD-EXPANSION 
SURFACE  EERJ:  VARIATIONS  OF  GROSS  THRUST  AND  NOZZLE 
LIFT  COEFFICIENTS  (BASED  ON  CAPTURE  AREA  A 0)  AS 
NOZZLES  ARE  ROTATED  DOWNWARD 
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is  much  greater  than 

(A„  =  A_){cos  0  )  *  gage  correction  for  CRJ 

r  3  v  u 

Obviously,  if  the  ERJ  is  boattailed  as  sketched,  the  pressure  on 
the  boattail  will  be  less  than  P^,  and  it  could  be  made  still  smaller  by 
streamlining  as  suggested  by  the  dashed  line;  moreover,  the  cowl  will 
add  to  the  lift.  (Cowl  lift  was  included  in  Figs.  12  through  15. )  If  the 
conventional  ramjet  is  pod- mounted,  then  the  pressure  below  its  nozzle 
v/ill  also  be  greater  than  the  pressure  above  its  nozzle,  so  that  the  ex¬ 
ternal  surfaces  will  increase  the  overall  lift  and  reduce  the  net  thrust. 

It  is  likely  that  this  effect  would  penalize  a  vehicle  employing  pod- 
mounted,  three-dimensional  conventional  engines  v/ith  vectored  nozzles, 
whereas  the  cowl  serves  to  increase  the  I./D  of  the  two-dimensional  1  RJ 

at  low  to  moderate  Pc  /?„.  If,  on  the  other  hand,  the  nozzle  is  ins'  e  a 

o  U 

fuselage,  these  pressure  effects  probably  w,  uld  be  eliminated.  Thus 
the  relative  merits  of  the  two  vectored -nozzle  engines  with  respect  to  lift 
and  thrust  will  depend  on  how  the  engines  are  integrated  with  airframes. 
Further  discussion  on  the  integration  of  the  EERJ  with  airframes  is 
given  in  a  later  section. 

Substitution  of  Simple  Wedge  Effusor 

In  connection  with  the  foregoing  calculi  ions  with  ideal,nonreflect- 
ing  expansion  surfaces,  it  was  pointed  out  that  for  Pq>  10,  corre¬ 
sponding  to  flow  turning  angles  less  than  one  radian,  Ag /A^  becomes 
small  {<  1  /2),  so  that  the  rear  portion  of  the  cowl  becomes  large  even 
though  the  front  was  of  the  minimum  dimensions.  Such  a  cowl  con¬ 
tributes  negatively  to  lift  while  still  introducing  an  appreciable  drag. 
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In  order  to  permit  operation  at  smaller  turning  angles  with  a  small,  low- 
drag,  low-cooling-requirerrent  cowl,  it  might  be  desirable  for  some  mis¬ 
sions  to  compromise  the  thrust  by  substituting  a  simple  wedge  in  place 
of  the  iueal  effusor.  A  wedge-effusor  engine  would  be  simpler  to  con¬ 
struct,  more  adaptable  to  geometry  variation  to  vary  B^  of  A^,  and 
would  provide  more  packaging  volume  under  the  effusor  and  make  it 
easier  to  control  the  c.g.  to  overcome  pitch. 

Figure  17  illustrates  the  flow  model  for  wedge  expansion.  This 
sketch  is  oversimplified,  because  the  Mach  lines  originating  from  the 
Prandtl-Meyer  expansion  around  the  cowl  lip  are  actually  curved  before 
they  reach  the  ramp,  as  a  result  of  refraction  by  earlier  Mach  lines 
which  have  already  been  reflected  from  the  surface.  This  may  be  seen 
in  Fig.  18,  which  shows  a  portion  of  the  graphical  characteristics 
solution  (obtained  by  the  method  described  in  Ref.  25)  for  each  of  three 
cases. 

By  comparing  various  cases  in  Figs.  18  and  19  one  can  see  two 
effects  which  further  enhance  the  refraction  effect.  Comparison  of  the 
ramp  intercepts  in  Fig.  18  for  the  solid-line  and  dashed-  line  character¬ 
istics,  representing  incorrect  use  of  the  frozen  sound  speed  a^  and  cor¬ 
rect  use  of  the  equilibrium  sound  speed  ae,  respectively,  shows  that 
the  latter  characteristics  have  been  refracted  more  for  a  given  degree 
of  turning  at  the  cowl  lip,  and  that  the  difference  increases  with  distance. 
The  second  effect  is  due  to  greater  refraction  as  equivalence  ratio  is 
increased,  as  may  be  seen  by  comparing  the  solid-line  intercepts  for 
E.  R.  =  1.  0  to  the  dotted-line  intercepts  for  E.R.  =  0.5.  The  result 

on  the  normal  force  coefficient;  F_  /q_A,_,  where  F.  is  the  total 

n  n 

normal  force  owing  to  expansion  over  a  ramp  of  length  i  is  shown  for 
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Fig  .  18  REAL  GAS  CHARACTERISTICS  FOR  EXPANSION  OVER  SIMPLE  WEDGE, 
ILLUSTRATING  REFRACTION  OF  MACH  LINES  AND  EFFECTS  OF  EQUIVALENCE 
RATIO  AND  OF  af  VERSUS  a  ON  REFRACTION 
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Fig.  19  COEFFICIENT  OF  NORMAL  FORCE  DOE  TO  EXPANSION  VERSOS 
EXPANSION  RAMP  LENGTH  FOR  SIMPLE  WEDGE — 

EFFECTS  OF  REFRACTION 
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these  three  cases  in  Fig.  19  for  the  range  of  ramp  lengths  of  interest. 

For  E.R.  =  1.  0  the  ambient  pressure  is  reached  at  i  =  1.52  h^ 

(where  h^  is  capture  height  corresponding  to  A^),  so  that  greater  lengths 
are  of  no  interest  for  the  equivalence  ratio.  Since  the  refraction  ^ect 
is  smaller  at  E.  R.  =  0.  5,  the  greatest  i/h^  of  interest  is  of  the  order 
of  5.  A  curve  for  E.  R.  =  0.  5  based  on  the  assumption  that  refraction 
can  be  neglected  completely  is  also  included  for  comparison  to  show  why 
the  wedge  results  obtained  in  this  manner  can  be  very  misleading. 

In  order  to  present  thrust  and  lift  forces  in  terms  of  A  / An,  the 

e  v 

following  relationship  is  employed: 


•? 

f 


? 

-  » 


i 

* 


5  I 


I 

I 

1 

I 

I 

1 


Ae/AQ  =  i  sin  94  +  (A4/AQ)  cos  04  (20) 

A  stream  thrust  analysis  shows  that  the  gross  thrust  will  now  be 
given  by 


F  =  F5  sin  e4+^  cos  e4 -&0  -  P0(Ae  -  AQ) 
n 


(21) 


where 


and 


F"  ~  C5  q0A0 


n 


n 


=  P4A4  (1  +  aft  2/R4T4) 


_  /  *• A  A  • 

e .  4  4 

4 


and  the  lift  force  on  the  configuration  will  be 


where 


L  =  L_  +  L.  +  L 
5  4  c 

net  net 


L5  ’  (F5  '  Po<)  003  94 
net  n 
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and  L  is  again  calculated  as  in  Appendix  I. 

Cnet 

The  results  are  summarized  in  Table  III.  The  highest  overall 

I.  L/D  calculated  was  5190  for  0=0.  524  rad  (30°)  in  the  series  in 
f  4 

which  was  held  constant  at  the  minimum  value  (determined  by  the 
F 

10°  diffuser  lip)  of  1.  008  A^.  For  this  turning  angle  of  30°,  P^/P^ 
is  near  unity  (approximately  0.  9),  so  that  no  further  expansion  would 
be  realistic.  The  value  of  5190  for  1^  L/D  is  only  75  percent  of  the 
maximum  value  cf  6900  calculated  for  the  isentropic  ramp  model  at 
this  equivalence  ratio. 

It  is  therefore  concluded  that  a  simple  wedge  effusor  would  not 
be  chosen  for  a  maximum  range  vehicle,  although  the  other  advantages 
mentioned  previously  might  well  lead  to  its  choice  for  short  range  mis¬ 
siles  or  for  higher  velocity  missiles  in  which  the  structural  and/or 
cooling  problems  associated  with  a  long  cowl  might  outweigh  the  loss 
in  propulsion  performance. 


4180 


X 


THC  JOHNS  HOPKINS  UNIWftfttTV 

APPLIED  PHYSICS  LABCHATORY 
SU.VK*  SPSINS  MA4YIAKO 


IV.  RANGE-WEIGHT  COMPARISON  OF  ERJ  WITH  OTHER  MISSILES 

1500-POUND  USEFUL  PAYLOADS 

Figure  20  compares  predicted  weight-range  performances  for  four 
ballistic  rocket  alternates,  one  glide  rocket,  and  three  ramjet  alter¬ 
nates,  the  last  of  which  would  incorporate  external  burning.  The  assump¬ 
tions  used  to  calculate  these  curves  are  discussed  below. 


General  Assumptions 

In  assigning  useful  payload  weights  for  the  rockets,  it  was 
assumed  that  half  of  the  nose  cone  weight  is  payload  (war¬ 
head)  and  half  metal  parts  plus  protection  (see  also  item  7). 
For  the  ramjet,  because  of  scaling  requirements,  it  was 
necessary  to  restrict  the  analysis  to  one  useful  payload  (or 
warhead)  weight  of  1500  pounds.  The  weight  of  the  guidance 
for  the  ramjet  was  assumed  to  be  independent  of  range. 

Thus  the  ramjet  curves  would  vary  if  useful  payload  weight 
were  changed.  The  ramjet  weight  advantage  would  in¬ 
crease  for  bigger  payloads  or  decrease  for  smaller  pay¬ 
loads. 

All  rockets  follow  the  optimum  trajectory  for  maximum 
range,  and  all  ramjets  follow  the  Breguet  trajectory  for 
maximum  range  after  rocket  boost  to  Mach  5  at  100,  000 
feet. 

Ratio  of  initial  weight  to  burnout  weight  for  each  rocket 

stage  in  a  given  missile  is  constant,  but  I  for  the  second 

sp 
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(and  third)  stage  is  assumed  to  be  20  percent  greater  than 

the  effective  I  of  the  first  stage  (to  allow  for  greater 

thrust  at  higher  altitude),  so  that  the  velocity  increment 

for  the  second  (and  third)  stage  is  20  percent  greater 

than  the  velocity  increment  for  the  first  stage.  A  3  percent 

drag  allowance  is  made  for  the  first  stage  so  that  I  = 

SPj 

0.  97  I  ,  and  I  =  I  =  1.  2  x  0.  97  I  .To  achieve 
sp0  sp2  sp3  spo 

these  impulses,  it  is  assumed  that  the  first  stage  nozzle 

has  an  expansion  ratio  near  four  and  the  second  (and  third) 

stage  has  an  expansion  ratio  of  20  to  30. 

4.  Each  missile  has  an  initial  acceleration  of  4g:s. 

5.  The  basic  rocket  stage  of  the  near  future  is  assumed  to  be 

a  solid  propellant  rocket  of  Ic  =  250  lb-sec /lb  and  struc- 

£po 

ture/fuel  weight  ratio  a  =  0.  07.  These  values  are  used 
for  curves  A,  B,  C,  and  E  and  for  the  boosters  for  ramjet 
curves  F,  G,  and  H. 

6.  The  probable  practical  limit  for  liquid  propellant  stages  of 

the  near  future  is  assumed  to  be  I  =  300  lb-sec /lb 

Spo 

(estimated  maximum  delivered  1  for  fluorine-hydrazine) 

sp 

and  a  =  0.  07.  This  a  value  is  believed  to  be  the  realistic 
goal  for  rockets  of  the  sire  indicated  by  Curve  D  (that  is, 
rockets  much  smaller  than  Atlas)  when  realistic  allowances 
are  made  for  the  minimum  control  weight,  fuel  residues, 
and  pressurization  residuals  which  must  be  assigned  to 
each  stage.  Indeed,  some  writers  now  feel  that  a  for  solid 
rockets  may  in  the  future  be  less  than  a  for  liquid  rockets 
of  comparable  total  impulse.  A  factor  contributing  to  this 
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belief  is  that  higher  solid  propellant  density  assures  higher 
thrust  per  unit  volume  for  the  solid  rocket. 

7.  For  the  glide  rockets  of  Curve  E,  and  L/D  of  4  was 
assumed  which  may  be  "about  right"  if  a  second  stage 
such  as  suggested  by  Eggers,  et  al. .  (see  Fig.  12,  Ref. 

26),  proves  to  be  practical.  This  suggested  design  is 
simply  a  conical-nosed  body  with  adjustable  petal-type 
control  surfaces  at  the  after  end.  The  vehicle  would  have 
relatively  low  drag  during  climb  and  low  aerodynamic  heat¬ 
ing  during  glide.  However,  the  addition  of  even  these  sim¬ 
plest  variable  control  surfaces  would,  because  of  the 
severity  of  aerodynamic  heating  on  the  control  surfaces,  add 
appreciably  to  the  glide  vehicle  weight.  The  glide  vehicle 
weight  is,  of  course,  the  payload  of  the  rocket  necessary  to 
accelerate  the  overall  missile.  A  20  percent  increase  in 
the  metal  parts  or  coolant  load  of  the  glider  would  increase 
the  launching  weight  7  percent.  This  fact  was  considered 

in  computing  Curve  E.  The  burn-out  velocity  versus  range 
relationship  used  to  obtain  the  curve  was  taken  from  Eggers, 
et  al. 

8.  The  ramjet  curves  are  based  on  a  relative  weight  breakdown 
for  an  advanced  ramjet  missile,  with  structures  assignable 
to  the  payload  scaled  up  to  a  1500  pound  payload  weight  and 
with  cooling  water,  dive  fuel,  and  related  structures  scaled 
to  total  ramjet  stage  weight.  One -half  pound  of  airframe 
weight  was  added  for  each  pound  of  cruise  fuel. 

9.  The  fuel  specific  impulses  chosen  for  the  ramjets  are  be¬ 
lieved  to  t .  realistic  for  kerosene  (I^.~  1500  lb-sec /lb  at 
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Mach  5)  and  HEF-3  (I^~  2000  lb-sec /lb  at  Mach  5).  These 
values  correspond  to  an  overall  thermal  efficiency  of  0. 5. 

On  the  assumption  of  a  nearly  constant  thermal  efficiency, 
the  product  of  VI^  remains  constant  for  higher  Mach 
numbers,  so  that  if  L/D  be  constant,  these  same  curves 
are  approximately  representative  for  higher  ramjet  cruise 
speeds,  for  example,  to  M  =  10. 

10.  The  lift-drag  ratios  assumed  for  the  ramjets  are:  5  for  a 
conventional  missile  at  Mach  5,  and  8  for  an  advanced  mis¬ 
sile  which  would  obtain  part  of  its  lift  from  external  burning. 
Since  an  1^  L/D  of  9000  seconds  was  estimated  for  external 
burning  of  kerosene  at  an  E.R.  of  0.  6  in  the  previous  section, 
the  present  use  of  L/D  =  2000  x  8  =  16,  000  seconds  for  a 
well -designed,  high-energy  fueled  engine  is  not  unreasonable, 
particularly  when  one  remembers  that  no  credit  for  engine 
weight  reduction  is  being  given  to  curve  H. 

Conclusion 


In  general,  the  conclusion  available  from  this  figure  is  that  for  a 
broad  and  important  span  cf  flight  distances  the  ramjet  has  an  appreciable 
gross-takeoff-weight  advantage  over  the  rocket  or  glide  rocket.  That  is, 
for  equal  payloads,  the  gross  takeoff  weight  of  the  ramjet  will  be  appre¬ 
ciably  less  than  that  of  either  the  glide  or  ballistic  rockets.  Possible 
comparisons  available  from  Fig.  20  are  shown  in  Table  IV. 
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-stage,  glide,  JP-5,  L/D  =  5  E  -  F  Ramjet  obtainable  with  current 

solid  technology  still  superior  at  ranges 

less  than  4000  n.  m.  Note  that 
gliding  descent  would  add  an  addi¬ 
tional  100  n.  m.  to  ramjet  range. 
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V.  DESIGN  PHILOSOPHY  FOR  HYPERSONIC 
EERJ  TRANSPORT  AIRCRAFT 

The  primary  advantage  which  an  external  expansion  ramjet  would 
have  over  a  conventional  ramjet  in  hypersonic  flight  is  the  heat  relief 
gained  through  radiation  from  the  surfaces  adjacent  to  the  burned  gases. 
In  addition,  a  signal  advantage  accrues  from  conversion  of  the  bursting 
forces  produced  in  a  ducted  nozzle  into  a  useful  component  of  lift  when 
the  burned-gas  expansion  takes  place  solely  on  the  under  side  of  the 
vehicle.  Consequently,  the  EERJ  must  be  an  integral  part  of  the  main 
body  or  wing,  and  no  pod-type,  thrust-generating  devices  should  be 
permitted  to  complicate  and  vitiate  the  basically  clean  design.  A  type 
of  flying  wing  is  thus  indicated. 


i 

\ 


I 

1 

I 


I 
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Apart  from  the  question  of  integrating  the  propulsion  and  aero¬ 
dynamic  features  of  the  hypersonic  ensemble,  the  greatest  obstacle  to 
practical  accomplishment  of  such  a  high-speed  vehicle  is  the  extreme 
range  of  flight  conditions  to  be  met.  It  is  assumed  that  the  transport 
must  take-off  and  land  on  its  own  power,  and  besides,  it  must  pass 
through  the  transonic  and  lower  supersonic  Mach  number  range  before 
the  variable  geometry  EERJ  can  become  operative  (at  about  Mach  3)  to 
provide  acceleration  to  approximately  Mach  7-10,  where  the  cruise 
flight  will  be  maintained. 

Since  there  is  neither  adequate  design  information  for  variable 
inlets  capable  of  operating  at  hypersonic  speeds  nor  any  backlog  of 
experience  on  combustion  processes  under  hypersonic  flight  conditions, 
no  unique  optimum  aircraft  design  can  be  visualized.  Three  versions 


■m  'mSM&6S*i»ZSH853^ 


-  87  - 


THC  JOHNS  HOSKINS  OHIVSMITY 

applied  physics  laboratory 

UVH  SWKNS.  MAffYlAMO 


of  hypersonic  jet  transports  ("idea  airplanes")  which  incorporate  turbojet 
engines  for  take-off  and  landing  and  ERJ  engines  for  high- speed  flight 
have  been  considered  qualitatively.  These  are  illustrated  in  Figs.  21 
through  24. 

Th'>  first  version  incorporates  a  delta  wing,  which  has  the  desirable 
property  of  exhibiting  very  little  center  of  pressure  movement  with  varia¬ 
tion  in  Mach  number  from  low  subsonic  to  '  ypersonic  speeds.  Inasmuch 
as  the  center  of  gravity  will  not  normally  change  its  position  markedly 
during  the  course  of  the  flight,  it  is  desirable  to  select  such  a  design, 
giving  constant  c.p.  location,  in  order  to  lessen  the  stability  and  con¬ 
trol  difficulties  which  would  otherwise  be  encountered.  However,  in 
or^er  to  achieve  low-speed  landing  this  delta- wing  version  depends  on 
large  wing  area,  lightly  loaded.  To  overcome  the  resulting  large  drag, 
the  turbojet  thrust  required  for  take-off  and  acceleration  will  be  large,  and 
required  turbojet  engines  become  dead  weight  and  volume  wasters  for  the 
high-speed  flight  portion  of  the  plane's  trajectory.  The  state  of  develop¬ 
ment  of  the  auxiliary  turbojets  thus  has  a  direct  bearing  on  the  ultimate 
design  of  a  truly  versatile  ERJ  hypersonic  design.  Figures  21  and  22 
show  take-off  and  cruise  configurations. 

Although  he  delta  wing,  in  the  absence  of  variant  combustion 
pressure  fields,  would  appear  to  be  an  admirable  choice  for  a  versatile 
multi-speed-range  aircraft,  it  may  be  found  that  the  stability  deteriorates 
under  burning  conditions.  Consequently,  the  other  line  of  approach  to 
the  attainment  of  trustworthy  dynamic  response  at  all  speeds  is  to  pro¬ 
vide  a  configuration  which  can  circumvent,  by  adequate  control,  the 
expected  severe  c.p.  changes.  The  second  version  has  a  tandem 
monoplane  (libellula)  configuration  which  has  the  advantage  of  equally 
disposed  lifts  (both  operating  in  the  effective  upward  direction)  which 


DELTA  PLANFORM  VEHICLE  IN  TAKE-OFF  CONFIGURATION 


-  72  - 


Fig.  24  COMPACT  "FLYING  WING"  VEHICLE  IN  TAKE-OFF  CONFIGURATION 
WITH  PECTORAL  WINGS  EXTENDED 
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may  be  adjusted  by  the  aid  of  spoilers  or  flaps  to  provide  a  very  large 
margin  of  total  airplane  c.p.  travel,  more  than  obtained  in  any  other 
design  having  a  feasible  tail  length.  The  slot  between  the  fore  and  after 
wings  also  acts  to  provide  ultra-high  circulation  for  slow  speed  landing. 
Severe  transition  problems  occurring  while  converting  from  low  to  high 
speed  .flight  configuration  may  be  anticipated.  Figure  23  is  an  upward 
view. 

The  third  model,  which  represents  attempts  to  eliminate 
fuselage  fairings  and  provide  a  more  compact  design  (passengers 
seated  within  the  wing  itself)  would  be  a  most  efficient  hyper  sonic -flight 
vehicle,  but  the  c.  p.  variations  over  the  speed  range  may  be  too  ex¬ 
cessive  to  be  dealt  with  by  conventional  controls.  Extra  wing  area  is 
provided  by  laterally  sliding  pectoral  flapped  planes,  and  this  fore-area 
control  versus  the  large  after  flap,  with  circulation-control  blowing- 
slots  at  its  entrant  edge,  may  be  sufficient  to  overcome  the  anticipated 
large  c.pw  changes.  Figure  24  is  a  3/4  front  view,  which  shows  the 
pectoral  flapped  planes  extended  for  landing. 

In  regard  to  the  requirement  for  thrust  generation  through  a  wide 
range  of  speeds,  it  is  believed  that  a  hypersonic  air-breathing  vehicle 
will  necessarily  incorporate  some  variable  geometry  at  the  expense  of 
structural  complexity.  The  delta-wing  model  illustrates  the  concept 
of  varying  the  shapes  and  positions  of  the  main  compression  and  expan¬ 
sion  surfaces  of  the  ERJ  while  holding  the  subsonic  diffusion  cowl  fixed. 
The  two  straight-wing  models  rely  on  the  supposition  that  the  air  flow 
area  between  the  compression-ramp-and-effusor  surface  and  the  out¬ 
rigger  vane  provides  an  adequate  control  variable,  which  is  varied 
simply  by  extending  the  vane  supports.  Although  this  expedient 
obviates  some  structural,  insulation  and  fuel- supply  problems,  it 
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implies  toleration  of  nonoptimum  jet  operation-  whether  off-design 
performance  of  the  EERJ  will  deteriorate  catastrophically  cannot  be 
predicted  at  present.  High  drag  in  the  base  section  of  the  EERJ  is 
also  to  be  expected  at  low  speed. 

In  any  event,  there  are  many  problems  associated  with  the  de¬ 
tails  of  the  flow  in  the  EERJ  and  with  such  full-gamut  speed- range 
aerodynamic  ensembles,  so  that  no  hard  and  fast  decisions  as  to  the 
main  features  of  such  a  hypersonic  vehicle  can  be  made  at  this  time. 
The  proposed  models  merely  suggest  that  there  are  apt  to  be  a  variety 
of  ways  in  which  the  ERJ  principle  can  be  wedded  to  hypersonic  aero¬ 
dynamic  shapes  to  provide  the  performance  capabilities  foreseen  from 
theoretical  anal yses. 
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APPENDIX  I 


Method  of  Estimating  Frontal  Area  A^,  Cowl  Wave  Drag 

Coefficient  C^  .Friction  Coefficient  C^, 
- D  - - - F 

c 

and  Lift  Coefficient  CT 


It  is  assumed  that  the  curved  parts  of  the  fore  and  after  cowl  sur¬ 
faces  approximate  cylindrical  surfaces  having  radii  equal  to  four  times 
the  local  duct  heights,  which  correspond  to  areas  A^  and  A^,  respectively. 
The  remaining  part  of  the  cowl  is  assumed  to  be  a  plane  surface  tangen¬ 
tially  extending  either  the  front  or  rear  cylindrical  surface,  depending  upon 
the  geometry  dictated  by  and  A^.  In  all  cases  it  is  assumed  that  the 
leading  edge  of  the  cowl  makes  an  angle  of  10°  with  the  flight  axis.  (See 
Fig.  25.) 
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Fig.  25  ASSUMPTIONS  ON  COWL  GEOMETRY 
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If  A^'  and  A^'  are  defined  as  the  minimum  projected  areas  required  to 
turn  the  flow  from  -10®  to  0®  and  from  0°  to  0 A,  respectively,  then 

A  '  =  4A2  (1  -  cos  (10®)) 

Ar'  =  4A4  (1  -  cos  e4) 

From  the  geometry  it  may  be  seen  that  when 

A  <  (1  +  A  '  -  A  A.  =  A  •  and  A  =  1  +  A  •  -  A 
e  fr  ft  r  fe 

but  when 

A  >  (1  +  AJ  -  A  '),  A  =  (A  +  A  1  -  1)  and  A  =  A  ' 
e  f  r  f  er  r  r 

The  maximum  cross-sectional  area  for  the  body  plus  cowl  is  the  frontal 
area: 

A_  =  1  +  A.  =  A  +  A 
F  f  e  r 

The  pressure  P  on  the  front  surface  or  lip  of  the  cowl  is  that  behind 
the  oblique  shock  for  =  5  and  the  10°  leading  edge  angle;  P^/Pq  = 
3.013. 

It  is  assumed  that  the  flow  separates  on  tne  rear  surface  and 

that  the  pressure  there  is  0.  5  P^.  With  these  assumptions,  the  cowl 

wave  force  F  is  estimated  by: 
c 

F  =  (P  -  F_)  A  -  (0.5  P  -  P_)  A 
c  cl  O  f  O  O  r 

=  (2.013  A.+  0.5A  )  P _ 
i  r  O 
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and  the  cowl  wave  drag  coefficient  is  given  by 


CD  =  Fc%\  -  '2-013Af+0.5Ar)Po/qoA 
c 

The  external  friction  drag  was  crudely  estimated  using  the  relationship 
for  laminar  boundary  layer  flow 


CfVv-  !•  2  at  Mach  5 

from  Hoerner,  Ref.  27  (p.  17-3).  The  Reynolds  number  per  foot, 

4 

R^/i,  was  calculated  to  be  1.5  x  10  for  the  free  stream  conditions  of 
395°R  and  0.  01  atmosphere  at  Mach  5;  it  was  assumed  that  to  a  first 
approximation  the  above  relationship  could  be  applied  directly  to  both 
the  top  surface  (R^  =  1.  5  x  104  *tQtal)  and  the  cowl  (R^  =  1.  5  x  104  *  j). 

and  that  the  sum  of  the  two  C^'s  so  calculated  would  approximate  for 
the  engine;  each  of  the  Cf's  was  related  to  engine  frontal  area  by  multi¬ 
plying  by  the  appropriate  length,  *total  or  *cowl>  divided  by  Ap. 

The  cowl  lift  force  is  estimated  from 


LC  =  <Pcl  “  P0)  [4A2  Sin  (10°)+  (Af  "  AfM  COt  (10°)] 
net 


+  (0.  5  P  -  P  )  [4A  sin  0  +  (A  -  A  ')  cot  0  .] 
U  O  4  4  r  r  4 


~  Lc  /qOAF 
net 
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APPENDIX  II 

Validity  of  Prandtl-Meyer  Expansion 
Calculations  of  Section  III 


In  Section  III  it  is  assumed  that  the  burned  gas  may  be  turned 
through  a  Prandtl-Meyer  expansion  to  the  point  0,.  =  0  (flow  parallel 
to  top  surface  of  missile)  for  final  pressures  as  low  as  ambient, 

Pc  /P  =  1.  At  the  same  time,  for  the  purpose  of  engine  drag  esti- 
mates,  it  is  assumed  that  there  is  a  separation  region  on  the  rear 
end  of  the  cowl  in  which  the  pressure  is  approximately  one-half  of 
the  ambient  pressure  P^.  In  order  to  be  satisfied  that  the  Prandtl- 
Meyer  expansion  calculations  are  valid,  the  following  models  have 
been  considered  for  the  interaction  between  the  burned  gas  and  the  ex¬ 
ternal  air  stream. 

In  all  probability  there  is  a  separation  region  or  recirculation 
zone  initiated  from  some  point  C  on  the  cowl,  as  shown  in  Fig.  26(a). 
The  first  case  to  be  considered  is  that  in  which  the  pressure  in  the 
recirculation  zone  decreases  from  Pq  at  point  C  to  some  lower  but 
constant  pressure  along  the  line  AB.  The  burned  gas  must  be  turned 
until  pressure  equals  the  pressure  along  AB,  and  the  external  air  flow 
must  also  be  turned  and  accelerated  so  that  its  pressure  matches  the 
decreasing  pressure  along  CB.  The  external  air  flow  and  the  burned 
gas  flow  meet  at  B;  both  flows  must  be  turned  to  equalize  pressures 
and  directions.  The  slip  line  between  the  flows  is  labeled  BD.  One 
point  to  be  clarified  is  whether  the  weak  shock  wave  BB'  required  to 
turn  the  burned  gas  flow  could  interfere  with  the  Mach  line  AA'  and 
thus  affect  the  calculated  expansion  force  on  the  effusor  surface. 
(Incidentally,  if  it  did,  a  considerable  increase  in  pressure  on  the 
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Fig.  26(a)  MODEL  USED  TO  EXAMINE  VALIDITY  OF  PRANDTL-MEYER  EXPANSION 
CALCULATIONS:  EFFECTS  OF  *1.0*  SEPARATlSrBESAli^  ” 
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effusor  surface  would  result,  compounded  to  some  degree  by  a  flow 

separation  in  the  effusor  boundary  layer  (see,  for  example.  Ref.  28), 

so  that  calculations  could  only  be  conservative  in  this  respec., )  This 

possibility  was  investigated  for  P_  / P._  =  5  and  P_  /P  =  1,  wht*  e  P 

5  0  5  0  5 

is  the  pressure  at  Mach  line  AA1,  and  the  burned  gas  flow  is  para,  s-i 
to  the  top  surface  of  the  engine  as  it  crosses  AA*.  These  pressure 
ratios  bracket  the  region  of  greatest  interest  with  respect  to  pro¬ 
viding  high  1^  L/'D.  In  both  cases  it  was  assumed  that  the  pressure 
along  AB  was  0.  5  P^  and  that  the  external  airflow  was  at  Mach  5  and 
Pq  at  point  C,  These  assumptions  required  the  external  air  to  be 
turned  through  5.  3°  (accelerated  to  Mach  5.  62)  upon  reaching  point  B 
and  then  deflected  16.  6°  (or  11.3°  below  the  flight  axis)  by  the  shock 
BB".  The  shock  angle  for  the  external  air  flow  at  Mach  5.  62  was  25°, 

and  the  pressure  along  BD  wafs  3.  23  Prt.  For  both  P_  /P^  =  5  and 

U  5  0 

P&/Pq  =  1  the  shock  line  BP1  and  the  Mach  line  AA!  were  divergent, 
so  that  there  could  be  no  interference. 


If  the  pressure  along  AB,  PAT1,  were  greater  than  P_,  then  a 

AJtjr  O 

shock  would  occur  at  point  C,  and  the  external  air  flow  again  would  be 
turned  and  compressed,  and  the  part  of  the  cowl  above  the  separation 


region  would  contribute  positively  to  engine  lift.  If  P_  were  also  a 
PAB’  then  BB"  and  wou^d  stil*  diverge.  For  the  unique  case 
P5  =  PAB  r  PQ  and  95  =  OV  the  line  CB  would  be  parallel  to  the  flight 
axis,  the  separation  region  and  slip  line  would  become  a  wake  at 


pressure  P^,  and  there  would  be  no  lift  on  the  rear  part  jf  the  cowl 
or  the  effusor. 


It  is  therefore  concluded  that  the  present  zero-angle-of-attack 
©5  =  O  calculations  are  valid  so  long  as  a  PAB*  At  positive  angles 
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of  attack,  the  lift  for'-*-  would  be  increased  (provided  that  the  diffuser 
efficiency  did  not  d  p  much),  and  the  thrust  force  decreased. 

The  extremely  unlikely  case  of  attached  external  flow  was  also 
considered.  (For  low  exit  flow  angle  04  or  for  positive  angles  of  attack, 
flow 'separation  could  be  delayed  or  prevented  by  suction  on  the  cowl 
trailing  surface,  but  in  the  absence  of  suction  or  for  flight  at  an  angle 
of  attack  less  than  8^.  it  is  quite  likely  that  separation  will  occur. )  If 
separation  does  not  occur,  then  in  these  zero-angle-of-attack  calcula¬ 
tions  the  external  air  flow  is  accelerated  U  a  very  high  Mach  number 
by  Prandtl-Meyer  expansion  around  the  cowl  and  then  deflected  by  a 
shock  wave  AB"  to  follow  the  slip  plane  AD  [Fig.  26(b)],  Approximate 
calculations  were  attempted  for  several  84's,  but  no  solution  was  found 
for  a  weak  shock  at  AB".  A  strong  shock  solution  is  also  impossible 
for  the  94's  of  interest,  because  the  pressure  in  the  subsonic  flow 
region  behind  the  shock  would  be  <  0.  01  P^,  and  the  burned  ,as  flo 
could  not  be  expanded  to  such  a  low  pressure  before  meeting  the  same 
shock.  The  separated  flow  model  is  therefore  the  proper  one  to  use. 
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Fig.  28(b)  MODEL  USED  TO  EXAMINE  VALIDITY  OF  PRANDTL- 
MEYER  EXPANSION  CALCULATIONS;  ATTACHED 
FLOW  BENEATH  COWL 
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